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ABSTRACT 


Streamlined  radomes  will  always  introduce  depolarization  because  of  the  fact  that  any 
material  surface  viewed  at  an  oblique  angle  has  different  transmission  coefficients  for  fields 
polarized  parallel  and  perpendicular  to  the  plane  of  incidence,  a  phenomenon  known  as  divorce. 
The  cross  polarized  fields  are  similar-  in  nature  to  those  responsible  for  Condon  lobes  of 
parabolic  reflectors,  and  which  are  known  to  make  those  antennas  susceptible  to  cross¬ 
polarization  jamming.  The  relative  strength  of  those  lobes  is  directly  related  to  the  amount  of 
divorce,  and  is  quite  sensitive  to  the  phase  difference  of  parallel  and  perpendicular  transmission 
coefficients.  A  mere  20°  of  phase  divorce  has  a  comparable  effect  to  3  dB  of  amplitude 
divorce. 

This  paper  presents  a  detailed  analysis  for  general  radome  shapes,  based  on  geometrical 
optics  (ray  tracing).  Example  results  are  presented  for  conical  and  tangent  ogive  radome 
geometries  used  in  conjunction  with  a  circular-aperture  antenna  that  is  mechanically  scanned. 
The  results  show  that  when  the  antenna  is  scanned  off  axis,  the  cross-polarized  pattern  changes 
from  a  Condon  lobe  structure  to  a  difference  pattern  in  the  direction  orthogonal  to  the  axis  of 
scan,  e.g.:  under  azimuth  scanning  the  cross-polarized  sum  pattern  resembles  an  elevation 
difference  pattern,  while  the  cross-polarized  elevation  difference  resembles  a  sum  beam.  Hence, 
interference  in  the  cross  polarization  will  introduce  a  tracking  error  in  the  direction  orthogonal 
to  the  scan  plane. 


ii 


CONTENTS 


LIST  OF  FIGURES . iv 

1.  INTRODUCTION .  1 

2.  METHODOLOGY  .  2 

2.1.  Radome  Analysis  Techniques .  2 

2.2.  Overview  of  the  Analysis  Method  .  2 

2.1.  Conical  Radome  with  Borcsightcd  Antenna .  4 

2.3.1.  Geometry  Description .  4 

2.3.2.  Transmission  Calculation .  7 

2.3.3.  Example  Results  .  10 

3.  DETAILED  ANALYSIS  PROCEDURE . 15 

3. 1.  Antenna  Scanning  -  Coordinate  Rotation . 15 

3.2.  Ray  Tracing  -  Intercept  Points . 18 

3.2.1.  Cone  . 18 

3.2.2.  Tangent  Ogive  . 19 

3.3.  Ray  Fixed  Coordinates  -  Parallel/Perpendicular  Decomposition  . 20 

3.4.  Transmission  &  Radiation  Calculation . 21 

4.  CALCULATED  RESULTS . 23 

4.1.  Comparison  of  Conical  and  Ogival  Radomes  . 23 

4.2.  Effects  of  Antenna  Scanning . 23 

5.  CONCLUSIONS  AND  RECOMMENDATIONS . 31 

REFERENCES . 33 

APPENDIX  A:  CODE  POCUMENATION . 34 


iii 


LIST  OF  FIGURES 

1 .  Generic  Radome  Geometry  with  Scanned  Antenna  .  3 

2.  Plane  of  Incidence  (Ray  Fixed  Coordinate  System) .  5 

3.  Conical  Radome  with  Boresighted  Antenna:  (a)  Perspective  View  Showing  Plane 

of  Incidence;  (b)  Nose-on  View  Illustrating  Polarization  of  Rays  from  a  Slot 
Array  Antenna  .  6 

4.  Surface  Plot  of  Cross-Polarized  Aperture  Amplitude .  9 

5.  Contour  Plot  of  Cross-Polarized  Aperture  Amplitude .  9 

6.  Sum  and  Azimuth-Difference  Patterns  (No  Radome) .  11 

7.  Co-  and  Cross-Polarized  Patterns,  15“  Phase  Divorce  . 12 

8.  Co-  and  Cross-Polarized  Patterns,  30°  Phase  Divorce  . 12 

9.  Co-  and  Cross-Polarized  Difference  Patterns,  15°  Phase  Divorce  . 13 

10.  Co-  and  Cross-Polarized  Difference  Patterns,  30°  Phase  Divorce . 13 

1 1 .  Contours  of  Constant  Main  Beam  to  Condon  Lobe  Ratio  (dB) 

vs.  Amplitude  Divorce  (dB)  and  Phase  Divorce  (deg) . 14 

12.  Sequence  of  Three  Coordinate  Rotations  for  Antenna  Scan  Calculation . 16 

13.  Test  Case  Geometries  (Cone  and  Ogive)  . 24 

14.  Surface  Transmissivity  vs.  Angle  for  Test  Case . 24 

15.  Co-  and  Cross-Polarized  Patterns  for  On-Axis  Scan: 

(a)  Tangent  Ogive  Radome;  (b)  Conical  Radome . 25 

16.  Cross-Polarized  Aperture  Amplitude  Distribution 

at  30°  Azimuth  Scan,  Ogive  Radome . 26 

17.  Cross-Polarized  Aperture  Phase  at  30°  Azimuth  Scan,  Ogive  Radome . 26 

18.  Antenna  Patterns  with  Ogive  Radome  at  30°  Azimuth  Scan: 

(a)  Cross-Pol  Pattern  Cuts;  (b)  Comparison  of  Co-  and  Cross-Pol  . 28 

19.  Ogive  Radome  Patterns  at  60°  Azimuth  Scan:  (a)  Cross-Polarized 

Patterns;  (b)  Comparison  of  Co-  and  Cross-Polarized  Pattern . 29 


iv 


20.  Antenna  Patterns  with  Ogive  Radome  at  30°  Azimuth  Scan: 

(a)  Elevation  Difference  Channel;  (b)  Azimuth  Difference  Channel  . 30 

Al.  Example  Input  File  RADOME2.INP  . 37 

A2.  Example  Transmissivity  Input  File  . 38 

A3.  Example  Output  File  RADOME2.APR . 39 

A4.  Example  Pattern  Output  File  RADOME2.PAT . 40 

A5.  Example  Output  File  RADOME2.SNM  (Surface  Normals) . 42 

A6.  Example  Output  File  RADOME2.IPT  (Intercept  Points)  . 42 

A7.  Example  Log  File  RADOME2.LOG  . 43 

A8.  Flowchart  for  Program  RADOME2 . 44 


v 


Radomes  used  to  protect  antennas  from  the  outside  environment  or  to  streamline  the 
platform  on  which  the  antenna  is  installed  are  almost  always  designed  to  appear  electrically 
invisible  to  the  extent  possible.  The  dominant  concerns  in  radome  design  are  maximizing 
transmissivity  and  minimizing  boresight  error.  A  third  concern  that  has  not  received  as  much 
attention  is  radome-induced  depolarization.  The  antenna  beam  will  almost  always  pass  through 
the  radome  at  an  oblique  angle.  Whenever  it  does,  it  is  susceptible  to  depolarization  simply 
because  any  material  surface  viewed  obliquely  presents  different  transmission  coefficients  to 
fields  polarized  perpendicular  and  parallel  to  the  surface,  a  phenomenon  sometimes  referred  to 
as  "divorce’’. 

The  purpose  of  this  analysis  is  to  assess  the  degradation  in  antenna  performance  as  a 
function  of  radome  transmissivity.  The  derivation  presented  here,  based  on  geometrical  optics, 
is  valid  for  most  radome  shapes  and  antenna  scan  directions.  The  results  are  applicable  to  any 
radome  material,  including  multilayered  dielectrics,  as  long  as  the  transmissivity  vs.  angle  of 
the  flat  surface  is  known. 

This  paper  will  first  present  a  simplified  analysis  for  a  conical  radome  with  an  on-axis 
antenna  beam,  in  order  to  illustrate  some  of  the  principles  involved.  Second,  the  complete 
analysis  for  arbitrary  radome  shapes  and  arbitrary  scan  directions  is  given,  with  illustrating 
examples  for  a  tangent  ogive  radome.  These  analyses  show  that  in  general,  a  radome.  that  has 
a  significant  level  of  divorce  will  generate  cross  polarized  patterns  that  are  remarkably  similar 
to  the  "Condon  lobes"  observed  with  parabolic  reflectors. 
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2.1.  Radome  Analysis  Techniques 

There  are  two  widely-accepted  numerical  methods  for  radome  transmission/  reflection 
analysis.  One  method  calculates  the  field  distribution  due  to  the  antenna  on  the  inner  radome 
surface  at  a  number  of  discrete  points,  and  the  far  field  as  an  integration  over  the  outer  surface, 
accounting  for  transmissivity  effects  [1],[2].  This  is  generally  referred  to  as  "surface 
integration."  It  is  the  most  accurate,  but  also  the  most  time-consuming.  It  is  the  method  of 
choice  when  sidelobe  levels  must  be  known  accurately.  The  second  method  is  ray  tracing, 
where  rays  emanating  from  the  antenna  are  traced  through  the  radome  to  an  equivalent  aperture 
whose  field  distribution  is  used  as  a  secondary  source  for  the  far  field  [3].  Ray  tracing  is 
computationally  faster  than  surface  integration.  Although  it  does  not  accurately  predict  sidelobe 
levels  outside  those  nearest  the  main  beam,  it  nonetheless  accurately  characterizes  the  properties 
of  the  main  beam,  and  those  of  the  difference  channel  beams  used  for  monopulse  tracking. 
Since  these  are  the  parameteis  of  greatest  interest  here,  the  ray-tracing  method  was  the  one 
chosen  for  this  study.  Moreover,  many  of  the  geometrical  derivations  used  in  ray  tracing  will 
also  be  required  for  the  surface  integration,  so  many  of  the  derivations  that  follow  would  be 
useful  in  developing  a  more  general  radome  simulation  based  on  surface  integration. 

2.2.  Overview  of  the  Analysis  Method 

Figure  1  illustrates  an  antenna  scanned  to  an  arbitrary  angle  inside  a  conical  radome. 
The  antenna  aperture  projects  through  the  radome  onto  an  equivalent  aperture  outside  the 
radome.  The  resulting  equivalent  aperture  has  an  amplitude  and  phase  distribution  that  are 
distorted  due  to  the  radome’s  imperfect  transmissivity  properties.  The  far  field  pattern  of  that 
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Figure  1.  Generic  Radome  Geometry  with  Scanned  Antenna 

distorted  aperture  distribution  gives  the  equivalent  properties  of  the  antenna  plus  the  radome. 
Thus,  the  problem  can  be  viewed  as  one  of  calculating  the  equivalent  aperture’s  amplitude  and 
phase.  Another  viewpoint  that  may  be  more  natural  from  an  optics  perspective  is  to  regard  the 
illuminated  portion  of  the  radome  as  a  "lens"  having  a  transmissivity  distribution.  Removing 
the  radome  and  multiplying  the  original  antenna  aperture  distribution  by  that  transmissivity  gives 
an  "equivalent"  antenna  that  includes  the  radome  effects. 

The  procedure  begins  by  establishing  a  series  of  sample  points  on  the  antenna  surface. 
If  the  antenna  is  an  array,  those  sample  points  are  simply  the  array  elements.  For  each  sample 
point,  there  is  a  ray  that  passes  through  the  radome  and  terminates  at  a  similar  point  in  the 
equivalent  aperture.  We  must  first  calculate  the  point  where  that  ray  passes  through  the  radome, 
and  the  direction  of  the  radome’s  inward  surface  normal  at  that  intercept  point.  Then  the 
incident  field  can  be  decomposed  into  components  parallel  and  perpendicular  to  the  plane  of 


incidence,  which  are  acted  upon  separately  by  the  two  components  of  the  transmissivity.  Then 
the  transmitted  field  components  are  recombined  into  components  co-polarized  and  cross- 
polarized  to  the  incident  field.  The  details  of  this  analysis  method  are  given  in  section  3,  but 
first,  a  simplified  case  is  discussed  to  show  the  general  concepts. 

2.1.  Conical  Radome  with  Borcsighted  Antenna 

To  illustrate  the  effects  of  radome-induced  de-polarization,  we  begin  with  the  fairly 
generic  problem  of  a  conical  radome  over  a  circular-aperture  antenna  that  looks  out  through  the 
cone’s  apex.  The  actual  antenna  properties  are  not  important  as  ’onp  as  we  can  assume  that 
antenna-radome  interactions  are  negligible.  Assuming  also  that  is  in  the  antenna’s 

far  field  then  all  rays  from  the  antenna  aperture  to  the  radome  are  _  »rai  1.  This  allows  us  to 
treat  the  radome’s  properties  as  an  alteration  of  the  antenna  aperture  distribution.  In  other 
words,  we  will  construct  an  equivalent  aperture  in  front  of  the  radome  onto  which  we  will 
project  the  antenna's  aperture  distribution  multiplied  by  the  radome’s  transmittance  function. 
This  equivalent  aperture  is  then  used  as  a  secondary  source  of  radiation  to  predict  the  far  field 
pattern  properties  of  the  antenna/ radome  combination. 

2.3.1.  Geometry  Description 

The  transmittance  function  depends  on  the  radome  surface  properties  and  on  the  incident 
wave  polarization.  For  instance,  the  incident  wave  may  have  its  electric  field  parallel  to  the 
surface,  referred  to  as  perpendicular  polarization  because  E  is  perpendicular  to  the  "plane  of 
incidence".  The  plane  of  incidence,  illustrated  in  Figure  2,  is  that  plane  containing  the  incident 
ray  and  'he  surface  normal  vector.  Perpendicular  polarization  is  sometimes  referred  to  as  TF7 
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Figure  2.  Plane  of  Incidence  (Ray  Fixed  Coordinate  System) 


polarization  when  discussing  flat  surfaces  whose  z  unit  vector  is  in  the  normal  direction. 
Alternately,  the  incident  wave  may  be  parallel  polarized  (TMj),  with  H  parallel  to  the  surface 
and  E  parallel  to  the  plane  of  incidence.  Most  often  the  wave  must  be  wiitten  as  a  linear 
combination  of  the  two  polarizations. 

Suppose  that  as  in  Figure  3a,  the  antenna,  located  in  the  x-y  plane,  is  boresighted  along 
the  z  axis,  looking  straight  out  the  nose  of  the  cone.  Any  incident  ray  is  parallel  to  the  z  axis. 
The  plane  of  incidence  is  always  the  constant-^  plane  containing  the  ray’s  source  location  (in 

A  A 

p,<t>  coordinates).  The  unit  vectors  uj  and  ux  are  aligned  with  the  electric  Field  components  in 
the  p  and  </>  directions,  respectively. 

Figure  3b  shows  the  same  situation  view-ed  looking  down  onto  the  antenna  along  the 
radome  axis,  but  also  representing  the  antenna  as  an  array  of  slots  polarized  in  the  y  direction. 
Consider  a  ray  emanating  from  a  slot  on  the  -Fy  axis:  Its  polarization  is  entirely  parallel  to  its 
plane  of  incidence  (the  y-z)  plane.  On  the  other  hand,  a  ray  from  a  slot  on  the  +  x  axis  is 
polarized  perpendicular  to  its  plane  of  incidence  (the  x-z)  plane.  If  the  radome  has  phase 
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Figure  3.  Conical  Radome  with  Borcsighted  Antenna:  (a)  Perspective 
View  Showing  Plane  of  Incidence;  (b)  Nose-on  View  Illustrating 
Polarization  of  Rays  from  a  Slot  Array  Antenna 
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divorce,  then  one  ray  will  be  delayed  more  than  the  other  as  they  pass  through  the  radome,  but 
neither  one  will  experience  depolarization.  Now  consider  a  ray  originating  along  the  <£=45° 
plane:  Its  polarization  is  evenly  split  between  parallel  and  perpendicular.  Amplitude  divorce 
will  cause  one  component  to  be  attenuated  more  than  the  other,  and  the  emerging  ray’s 
polarization  will  be  different  than  that  of  the  incident  ray  (but  still  linear  polarization).  Phase 
divorce  will  cause  one  component  to  be  delayed  relative  to  the  other,  and  the  emerging  ray  will 
have  elliptical  polarization.  This  effect  is  the  cause  of  radome  depolarization. 


2.3.2.  Transmission  Calculation 

Knowing  the  angle  between  the  incident  ray  and  the  radome,  we  can  find  a  transmissivity 
for  each  of  the  two  polarizations.  Let  -y  y  (a)  and  y±(a)  denote  the  complex  transmissivity  for 
the  parallel  and  perpendicular  components  of  a  ray  incident  at  the  angle  a  from  the  radome’s 
inward  surface  normal.  For  this  special  case  of  conical  radome  with  boresighted  antenna,  every 
ray  from  the  antenna  has  the  same  incidence  angle  a.  Let  the  antenna  be  linearly  polarized  in 
the  y  direction.  Upon  reaching  the  radome’s  interior  surface,  the  incident  field  decomposes  into 
orthogonal  components: 


E[  ‘E‘y  u,  ; 


then  the  transmitted  components  arc 


E*  -  r< 


t  p  t 

»JL 


E;  =  E[ax-s  +  E'Lw 

E'X=E[UX-X  +  E'LU^X 


(2.1) 


(2.2) 


(2.3) 


Equations  (2. 1)  through  (2.3)  are  completely  general,  and  apply  to  any  radome  geometry  as  long 
as  the  incident  wave  is  y-polarized.  Now,  specializing  them  to  the  geometry  in  Figure  3: 


y  ■  tij  =  y  *  - p  =  -sin<f> 

A 

$  •  =  y  ■  ~4>  =  -cos<|> 

E\  =  'Ey  Y,sm4> 

E*  =  ~E*  yx  cos <J) 

Ey  =  E‘  ( y ,  sin2 <J)  +  yx  cos2<J>) 
Ex  =Ey  (Y|  -  YA)sm4>cos<i) 


(2.4) 


(2.5) 


(2.6) 


(2.7) 


From  the  last  expression  it  is  clear  that  a  cross-polarized  component  is  introduced  any 
time  there  is  a  difference  in  either  the  amplitude  or  the  phase  of  the  transmission  coefficients  (di¬ 
vorce).  That  cross  polarization  component  is  maximum  along  the  diagonals  tf>=450, 135°.  The 
surface  plot  of  Figure  4  shows  the  shape  of  this  cross-polarized  distribution  (the  original  antenna 
distribution  is  a  simple  cosine  taper).  This  is  the  amplitude  distribution  in  the  equivalent 
aperture  located  "in  front"  of  the  radome.  Figure  5  is  a  contour  plot  of  the  magnitude.  It  is 
extremely  important  to  note  that  each  quadrant  has  constant  phase,  but  alternating  quadrants  are 
J80°  out  of  phase.  Also  note  that  whenever  there  is  no  divorce,  there  is  no  cross  polarized 
component,  and  the  co-polarized  (y)  component  experiences  only  a  constant  phase  shift.  The 
form  of  the  cross  polarization  components  shown  above  are  identical  to  those  experienced  with 
parabolic  reflectors  and  known  to  be  the  source  of  " Condon  lobes "  that  made  those  antennas 
susceptible  to  cross-polarization  jamming  [4], (5: 153-155],  [6:337], 
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2.3.3.  Example  Results 

The  equivalent  aperture  distribution  is  now  treated  as  a  source  of  radiation,  as  though  it 
were  itself  an  antenna.  When  there  is  no  divorce,  the  far  field  patterns  in  the  diagonal  (<£=45°) 
plane  are  as  shown  in  Figure  6.  These  patterns  are  normalized  to  a  uniformly  illuminated 
aperture  with  unit  current  amplitude.  The  beam  peak  is  lower  than  0  dB  due  to  "taper  loss". 
The  sum  pattern  is  due  to  a  cosine  amplitude  distribution:  1  4-  cos(xp/a),  where  a  is  the 
aperture  radius.  The  azimuth  difference  pattern  is  formed  by  phase  shifting  the  left  side  of  the 
aperture  (x>0)  relative  to  the  right  side  by  180°,  and  setting  the  amplitude  in  the  center  of  the 
antenna  to  zero.  A  monopulse  tracking  antenna  will  usually  have  an  elevation  difference  as 
well,  formed  by  phase  shifting  top  and  bottom  halves  of  the  antenna  aperture.  The  difference 
lobes  have  equal  magnitude,  but  180°  phase  difference. 

The  form  of  the  cross  polarized  distribution  due  to  divorce  suggests  that  we  will  not  see 
any  cross  polarized  pattern  in  either  the  <£=0°  or  <£=90°  planes  because  of  its  anti-symmetry. 
For  instance,  measuring  along  the  horizontal  axis,  the  contributions  to  the  far  field  from  the  top 
and  bottom  of  the  aperture  are  equal  and  opposite.  The  most  pronounced  effects  are  in  the 
diagonal  plane,  where  those  contributions  tend  to  reinforce  each  other.  Thus,  measurements  of 
antenna  cross  polarization  must  be  taken  in  the  diagonal  scan  plane,  not  the  principal  planes 
when  the  antenna  is  boresighted  with  the  radome  axis. 

Figures  7  and  8  are  the  sum  pattern  in  the  diagonal  plane  when  the  radome’s  E  and  H 
plane  transmissivity  have  15°  and  30°  phase  divorce,  respectively,  and  no  amplitude  divorce. 
The  only  effect  on  the  co -polarized  pattern  is  a  slight  rise  in  the  sidelobes.  But  the  cross- 
polarized  term  is  significant  because  of  its  resemblance  to  the  difference  pattern.  It  actually  has 
four  main  lobes,  one  in  each  quadrant  with  nulls  in  the  principal  azimuth  and  elevation  planes. 
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ANGLE  (DEG) 

Figure  6.  Sum  and  Azimuth-Difference  Patterns  (No  Radome) 

Figures  9  and  10  are  the  cc-  and  cross-  polarized  difference  patterns,  also  with  15°  and 
30°  phase  divorce.  The  cross-pol  difference  pattern  is  very  similar  to  the  co-pol  and  its  main 
lobes  are  also  phase  shifted  relative  to  each  other.  The  effect  this  might  have  on  a  tracking 
system  depends  on  whether  those  lobes  have  the  same  phase  or  opposite  phase  of  the  co¬ 
polarized  difference  pattern  lobes. 

Figure  1 1  shows  contours  (in  dB)  of  the  main  beam  to  Condon  lobe  ratio  as  a  function 
of  both  amplitude  divorce  (abscissa)  and  phase  divorce  (ordinate).  It  illustrates  that  the  radome’s 
cross-pol  performance  is,  if  anything,  more  sensitive  to  phase  divorce  than  to  amplitude  divorce. 
For  example,  20°  of  phase  divorce  is  comparable  to  3  dB  of  amplitude  divorce!  (However, 
large  amplitude  divorce  also  implies  strong  internal  reflections,  which  may  pose  more  serious 
problems  than  depolarization.)  Figure  1 1  may  be  used  to  specify  maximum  tolerances  on  the 
radome  properties  in  order  to  permit  a  desired  cross-polarization  level. 
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Figure  7.  Co-  and  Cross-Polarized  Patterns,  15°  Phase  Divorce 
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Figure  8.  Co-  and  Cross-Polarized  Patterns,  30°  Phase  Divorce 
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jure  9.  Co-  and  Cross-Polarized  Difference  Patterns,  15°  Phase  Divorce 
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Figure  10,  Co-  and  Cross-Polarized  Difference  Patterns,  30°  Phase  Divorce 


Figure  11.  Contours  of  Constant  Main  Beam  to  Condon  Lobe  Ratio  (dB) 
vs.  Amplitude  Divorce  (dB)  and  Phase  Divorce  (deg) 

In  many  respects,  the  cone  shape  is  the  worst  for  cross  polarization  when  the  antenna  is 
not  scanned  because  the  incidence  angle  is  very  steep  where  the  strongest  part  of  the  beam 
passes  through.  The  boresight  scan  case  is  also  not  representative  of  a  typical  scenario  because 
the  depolarization  at  off-boresight  scan  angles  might  be  less  because  the  incidence  angle  is 
shallower.  Thus  it  is  clear  that  a  more  detailed  analysis  is  necessary  to  generalize  the  above 
results  to  more  practical  situations.  The  following  section  extends  this  analysis  method  to  a 
more  realistic  radome  shape  (tangent  ogive),  and  with  the  antenna  scanned  away  from  nose-on 
incidence. 
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3.  DETAILED  ANALYSIS  PROCEDURE 


This  section  expands  on  the  generic  analysis  procedure  used  for  the  above  conical  radome 
analysis.  Analyses  for  off-axis  scan  are  presented  for  both  conical  and  tangent  ogive  radomes. 
The  antenna  is  considered  to  be  a  circular  flat-plate  array  that  is  mechanically  scanned.  The  end 
product  from  this  analysis  is  a  general  simulation  program,  RADOME2,  described  in  the 
Appendix.  A  brief  outline  of  the  procedure  is  as  follows: 


(1)  Find  locations  of  rotated  antenna  aperture  sample  points  - 
this  establishes  the  source  and  destination  points  for  rays 
passing  through  the  radome; 

(2)  Find  the  points  where  the  rays  intersect  the  radome; 

(3)  Decompose  the  incident  field  into  polarization  components; 

(4)  Determine  transmission  coefficients; 

(5)  Combine  polarization  components  to  find  equivalent  aperture 
distribution;  and 

(6)  Calculate  the  radiated  fields. 


3.1.  Antenna  Scanning  •  Coordinate  Rotation 

The  antenna  wilt  be  assumed  to  be  an  array  whose  elements  are  distributed  over  a  circle 
of  radius  a  in  the  z=0  plane.  When  the  antenna  is  mechanically  scanned  to  the  angle  0o  from 
the  z  axis  in  the  4>0  plane  (polar  coordinate  angles),  most  elements  will  be  shifted  either  up  or 
down  in  z,  and  their  x  and  y  coordinates  will  also  change. 

In  order  to  find  the  rotated  element  locations  (xr,yr,zr)  in  terms  of  the  original  locations 
(x,y,z),  we  will  view  the  antenna’s  scanning  as  a  sequence  of  three  coordinate  transformations: 
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(a)  R, 


Figure  12.  Sequence  of  Three  Coordinate  Rotations  for  Antenna  Scan  Calculation 

(1)  Rotation  by  <t>0  about  the  z-axis,  producing  coordinates  (x’,y’,z’); 

(2)  Rotation  by  60  about  the  y’  axis,  resulting  in  (x”,y”,z”);  and 

(3)  Rotation  by  -<f>0  about  the  z”  axis,  resulting  in  (xr,yr,zr). 

These  coordinate  rotations  are  illustrated  in  Figure  12.  Each  of  these  transformations  may  be 
characterized  by  a  rotation  matrix  [7:199-200]: 


cos4»o  0 

cos0o  0  -sin0o 

cos«|>0  -sin4>0  0 

-sin4>0  cos4»0  0 

0  1  0 

sin<|>0  cos<J>0  0 

0  0  1 

sin0o  0  cos0o 

rH 

o 

o 

and  the  new  coordinates  are  found  through  matrix  multiplication: 
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(3.2) 
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>, 

=  RlR2Ri 

y 

z 

cos2  4>0  cos  0O  +  sin2  4>0 
cos4>0sin<))0(l  -cosGq) 
cos4>osin0o 


cos4>0sin4>0(l  -cosOq) 
sin2  $ocos0o  +  cos2  $0 
-sin<J>osin0o 


-cos<{>osin0o 
sin<|)osin0o 
cos  0Q 


(3.3) 


and  using  the  fact  that  z=0,  we  finally  obtain 


xr  =  x(cos24>ocos0o  +  sin2<t>o)  +  ycos<J>0sui<t>0(l -cosO^ 
yr  =  jccos4»0sin<i>0(l -cosGq)  +  y(sin2<j>ocos0o +  cos2<J>0) 
Zr  =  xcos<J>osin0o  -  ysin«J>osin0o 


Eqs.  (3.4)  give  the  coordinates  of  points  on  the  rotated  antenna.  The  ray  from  each  point  is 
described  by  the  line  joining  the  point  (x^y^z,)  to  a  corresponding  point  (xa,ya,za)  on  the 
equivalent  aperture: 

xa  =  xr  +  /isin60cos(t>0 

ya  =  yr+hsw%swi>0  (3.5) 

za  =  zr  +  h  cos  9q 

where  h  is  the  radome  height  measured  from  the  z=0  plane  to  the  tip.  Finally,  the  ray  is 
described  by  the  parametric  equation: 

r(0  -  r„  +  X(xa-xr)t  +  y(ya-y,)t  +  z(z0-zr)t  O^r^l  (3-6) 
The  next  step  is  to  find  the  point  where  this  ray  intersects  the  radome  surface.  The  following 
will  describe  this  calculation  for  a  cone  and  for  an  ogive. 
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3.2.  Ray  Tracing  -  Intercept  Points 

3.2.1.  Cone 

The  general  form  for  an  equation  describing  a  three-dimensional  surface  is  g(x,y,z)=0. 
For  a  cone  with  base  radius  b  and  height  h: 


gixyj)  =  x2+y2-b2(l-^)2  (3.7) 

h 

The  point  where  the  ray  (3.6)  intercepts  the  radome  is  found  by  substituting 

xr+t(xa-xf),  yr+t(ya-yr)  and  Zr+tfz^z,)  for  x,y  and  z,  respectively,  in  (3.7).  The  result  is  the 

following  quadratic  equation  g(t)-0: 

g(t )  =  At 2  +  Bt  +  C 

*  =  +  O'.-?,)2  -  ^u.-*,)2 

B  -  2[(jca -xr)xr  *  ( y„-yr)yr  *  ^r(z0-z,)(/i-z,)]  {i"j) 

hr 

C  =  x2  +  y)  -  -Ah-z,)2 
h 

After  solving  (3.8)  for  t,  the  coordinates  of  the  intercept  point  (xt,yt,zv)  are  found  by  substituting 
t  into  (3.6).  Note  that  the  quadratic  gives  two  solutions  for  t,  but  the  one  with  0<t^  1  is  the 
only  correct  choice. 

We  will  also  need  to  know  the  surface  normal  at  the  point  of  intersection.  The  surface 
normal  anywhere  on  the  surface  is  in  the  direction  of  the  gradient  of  g: 


Vg  =  *(2x)  *  y(2y) 


(3.9) 


The  normal  at  the  intercept  point  is  simply  (3.9)  evaluated  at  (x^y^Zj).  The  angle  of  incidence, 
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a,  is  found  from  cos(a)=(k*n)/|k|  where  k  is  the  propagation  vector(k  =  ra  -  rr)  and  n  is  the 
surface  normal. 


3.2.2.  Tangent  Ogive 

An  ogive  is  simply  a  surface  of  revolution  formed  by  two  intersecting  circular  arcs.  A 
tangent  ogive  is  a  special  case  in  which  the  base  is  perpendicular  to  the  sides  everywhere  they 
intersect.  Then  if  the  base  of  the  ogive  is  joined  to  a  cylinder  with  the  same  radius,  both 
surfaces  are  tangent  a*,  the  points  where  they  join.  A  general  equation  for  the  tangent  ogive  with 
base  b  and  height  h  is 

gixyj.)  =  x2  +  y2  -  ()/R2-z2  -of  =  0 

R  h2  +  b2  (3.10) 

2b 

a  =  R-b 

Substituting  the  components  of  the  parametric  vector  (3.6)  into  (3.10)  produces  the  following 
expression  g(t)=0  to  be  solved  for  t: 

g(t)  =  At2  +  Bt  +  C  +  2a\J  Dt2  +Et+F  =  0 

A  =  (Xa~Xr)2  +  (yo~yr)2  -  (Za~*r)2 

B  =  2[(xa-xr)xr  +  ( ya-yr)yr  +  ( ?a~zr)zr ) 

C  =  x)  +  y)  +  z]  -  (R2+a2)  0.H) 

D  -  za-zr )2 
E  -  -2zr(za-zr) 

F  -  Rl~z) 

Unfortunately,  this  is  not  a  simple  quadratic  as  in  the  case  of  the  cone,  so  it  must  be  solved 
iteratively  using  Newton’s  method: 
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g(ft) 

s' Ok) 


(3.12) 


*t.i  h 


The  final  converged  value  for  t  is  substituted  into  (3.6)  to  find  the  intercept  coordinates  (x^y,,^. 
The  surface  normal  is  the  gradient  evaluated  at  (Xj.yt.z^: 


Vg  =  x(lx)  *  y(2 y) 


+  i( 2z) 


V 


(3.13) 


3.3.  Ray  Fixed  Coordinates  -  Parallcl/Perpcndicular  Decomposition 

The  properties  of  the  radome  shell  material  arc  most  conveniently  characterized  in  terms 
of  parallel  and  perpendicular  transmissivities.  The  ray-fixed  coordinate  system,  illustrated 
earlier  in  Figure  2,  is  used  to  express  the  incident  and  transmitted  fields  in  terms  of  parallel  and 
perpendicular  components.  Referring  to  that  figure,  the  vector  k  is  in  the  direction  of 
propagation:  ic  =  ra  -  rr.  At  the  point  this  ray  intersects  the  radome,  the  radome’s  surface 
normal  is  n.  The  unit  vector  ux  is  perpendicular  to  both  k  and  n  since  both  of  the  latter  are  in 
the  plane  of  incidence.  It  is  given  by 

=  JJLL  (3.i4) 

then  the  orthogonal  unit  vector  u  j  is  found  from 

(3.15) 

The  three  unit  vectors  ux ,  U|  and  k  form  the  basis  vectors  for  the  ray  fixed  coordinate  system. 
The  incident  field  is  decomposed  using 


20 


(3.16) 


e{  =  Ei-ai 

E[  =  El  &x 

3.4.  Transmission  &  Radiation  Calculation 

The  two  components  of  incident  electric  field  parallel  and  perpendicular  to  the  plane  of 
incidence  are  scalars.  The  transmitted  field  components  are  simply  the  product  of  those  scalars 
and  the  corresponding  complex  transmissivities  yj  and  7X.  The  transmitted  field  is 


E‘  =  +  (E[yx)ux 


=  E[ux 


(3.14) 


This  transmitted  field  is  resolved  into  co-polarized  and  cross-polarized  components  as  follows: 
Let  e'  be  a  unit  vector  in  the  direction  of  the  incident  electric  field,  then  the  unit  vector  in  the 

A  •  A  •  A 

incident  magnetic  field  direction  is  h1  =  -c'  x  k.  The  transmitted  components  are 


E'cp  =  E'  -il 
E[p  =  E'  ■£' 


(3.18) 


These  scalar  components  are  regarded  as  the  co-  and  cross-polarized  excitations  of  the  equivalent 
aperture. 

The  antenna’s  co-polarized  radiation  pattern  is  calculated  from 


£cp(0,< M  =  f;£'f(m)g(0,4.)«v“'-™6cO!*,,,-“esM)  (3.19) 

m  =  l 

where  f (0,<£)  is  the  radiation  pattern  of  an  array  clement.  The  complex  coefficients  Elcp(m) 
include  the  amplitude  weighting  for  low  sidclobcs  as  well  as  beam  steering  phase  shifts,  if 
applicable.  The  subscript  m  denotes  the  index  of  the  antenna  element  where  the  ray  originated 
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and  (xm,ym)  are  its  coordinates  before  scanning.  In  other  words,  the  pattern  computation  is 
always  done  relative  to  the  antenna’s  boresight  direction,  not  the  radome’s.  The  cross-polarized 
radiated  field  calculation  is  the  same  as  (3.19)  with  subscripts  XP. 

This  completes  the  analysis  procedure  derivation.  The  program  RADOME2,  described 
in  the  Appendix,  incorporates  the  above  equations  (3.1)  through  (3.19).  It  is  capable  of 
arbitrary  scan  angles,  cone  and  tangent  ogive  radome  shapes  with  arbitrary  height  and  base 
radius  and  mechanically  or  electronically  scanned  antennas.  The  initial  antenna  amplitude 
distribution  is  a  Taylor  distribution  with  user -specified  sidelobe  level.  The  radome  trans¬ 
missivity  as  a  function  of  incidence  angle  is  entered  through  a  data  file,  and  the  program 
interpolates  to  find  the  transmissivity  at  the  required  incidence  angles.  The  following  section 
presents  some  of  the  results  from  the  program. 
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4.1.  Comparison  of  Conical  and  Ogival  Radomes 

The  test  case  geometry  is  shown  in  Figure  13.  The  radome  has  a  base  radius  of  7.5X 
and  a  height  of  30X.  The  antenna  has  an  aperture  radius  of  5X.  The  radome  shell  has  the 
transmissivity  vs.  angle  shown  in  Figure  14.  This  is  a  typical  characteristic  for  a  multilayer 
radome:  the  amplitude  is  well  matched  out  to  the  design  angle  (here  45°)  after  which  the 
components  begin  to  separate,  while  the  phase  begins  to  separate  immediately  away  from  normal 
incidence  with  the  difference  increasing  nearly  linearly  with  incidence  angle.  The  initial  antenna 
distribution  is  a  Taylor  distribution  (specialized  to  a  circular  aperture  [8])  with  n =4,  meaning 
that  the  first  4  sidleobes  on  each  side  of  the  main  beam  are  maintained  at  -40  dB  relative  to  the 
main  beam,  and  the  remaining  sidelobes  are  progressively  lower. 

It  is  clear  from  Figure  13  that  when  the  antenna  beam  is  directed  through  the  radome 
apex,  the  incidence  angles  arc  lower  for  the  ogive  shape  than  for  the  cone  shape.  The  result  is 
that  the  ogive  (Figure  15a)  has  much  lower  Condon  lobes  than  the  cone  (Figure  15b).  However, 
the  on-axis  scan  is  the  worst  case  for  the  cone,  and  the  incidence  angle  decreases  rapidly  as  the 
antenna  beam  scans,  so  that  in  the  range  of  30°-60°  scan,  its  performance  is  comparable  to  that 
of  the  ogive. 

4.2.  Effects  of  Antenna  Scanning 

The  final  set  of  calculations  involve  scanning  the  antenna  in  the  <£=0  plane,  that  is  in  the 
azimuth  plane  corresponding  to  y=0.  The  antenna  is  vertically  polarized,  with  electric  field  in 
the  y  direction.  At  30°  scan,  the  cross  polarized  aperture  amplitude  is  as  shown  in  Figure  16. 
It  is  now  symmetric  only  about  the  x  axis,  unlike  the  on  axis  case  (see  Fig.  5),  which  was 
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(deg) 


symmetric  about  both  axes.  The  cross-polarized  phase,  shown  in  Figure  17,  is  antisymmetric. 
The  result  of  a  symmetric  amplitude  and  antisymmetric  phase  is,  of  course,  a  difference  pattern. 
Thus,  there  are  only  two  lobes  in  the  cross-polarized  pattern  instead  of  four.  Those  lobes  appear 
in  the  elevation  plane  as  well  as  the  diagonal  plane,  as  shown  in  Figure  18a.  The  cross- 
polarized  azimuth  pattern  is  not  quite  zero  because  of  a  slight  asymmetry  in  the  element 
locations.  Figure  1 8b  shows  two  of  the  co-polarized  patterns  for  amplitude  comparison.  Figures 
19a  and  19b  are  the  same  patterns,  but  with  the  antenna  scanned  to  60®.  Scanning  in  the 
elevation  plane  has  a  similar  effect,  except  that  the  cross-polarized  patterns  will  have  twin 
difference  lobes  in  the  azimuth  pattern  instead  of  the  elevation  pattern. 

Last,  Figure  20  shows  the  response  of  the  two  difference  channels  for  the  same  case  as 
Figure  18,  i.e.  30°  azimuth  scan.  The  cross-polarized  pattern  in  the  elevation  difference  channel 
looks  more  like  a  sum  beam.  The  azimuth  difference  channel,  on  the  other  hand,  still  has  a  null 
in  the  antenna’s  broadside  direction. 

The  anticipated  effect  of  this  is  that  when  the  antenna  is  looking  out  through  the  left  or 
right  side,  the  jammer  may  cause  an  error  in  the  elevation  direction;  and  when  the  antenna  is 
looking  out  through  the  top  or  bottom  of  the  radome,  a  cross-polarized  jammer  may  cause  a 
tracking  error  in  the  azimuth  direction.  The  exact  nature  of  the  tracking  error  and  the  J/S  ratio 
required  to  defeat  the  monopulse  tracking  depends  on  several  other  factors,  such  as  the  receiver 
design  and  the  nature  of  the  jamming  signal.  Those  are  recommended  topics  for  future  study. 


26 


RELATIVE  POWER  <dB)  RELATIVE  POWER  (dB) 


RELATIVE  POWER  (dB)  RELATIVE  POWER  (dB) 


RELATIVE  POWER  (dB)  RELATIVE  POWER  (dB) 


5.  CONCLUSIONS  AND  RECOMMENDATIONS 

Every  radome  will  generate  some  depolarization  simply  because  the  antenna  beam  passes 
through  it  at  oblique  angles,  and  any  material  surface  viewed  obliquely  has  different  parallel  and 
perpendicular  transmissivities  (divorce).  This  report  has  shown  that  the  cross  polarized  fields 
generate  difference  patterns  due  to  their  antisymmetric  nature.  An  on  axis  beam  will  generate 
the  classic  Condon  lobe  pattern  observed  earlier  in  parabolic  reflectors.  When  the  antenna  beam 
is  scanned  off  axis,  so  that  it  looks  out  through  the  side  of  the  radome,  the  cross  polarized 
pattern  has  only  two  lobes,  whose  beamwidths  and  locations  are  identical  to  those  of  the  mono¬ 
pulse  difference  pattern’s  main  lobes,  but  their  response  is  seen  in  the  receiver’s  sum  channel. 
The  result  is  a  susceptibility  to  cross-polarization  jamming. 

The  relative  strength  of  the  cross-pol  lobes  was  shown  to  be  directly  related  to  the 
transmissivity  properties  of  the  radome  structure.  The  transmission  coefficients  for  parallel  and 
perpendicular  polarizations  are  usually  different  in  both  amplitude  and  phase.  It  was  shown  that 
the  phase  divorce  is  of  greater  concern  than  the  amplitude  divorce,  with  a  mere  10°  phase 
difference  causing  more  depolarization  than  1  dB  amplitude  difference. 

This  study  used  a  fairly  simple  geometrical-optics  analysis.  The  results  are  expected  to 
be  reasonably  accurate  for  angles  near  the  main  beam.  It  is  recommended  that  these  results  be 
compared  with  measured  data  to  validate  their  accuracy.  The  sidelobe  levels  shown  are  not 
expected  to  be  accurate  because  they  do  not  account  for  internal  radome  reflections.  If  good 
estimates  of  sidelobe  levels  are  required,  we  recommend  that  the  analysis  method  be  modified 
to  use  surface  integration. 

This  report  identifies  a  potential  problem.  The  extent  of  a  system’s  vulnerability  to 
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cross-polarization  jamming  depends  on  the  receiver  design.  We  recommend  that  the  results 
presented  here  be  used  to  simulate  the  tracking  performance  of  a  monopulse  radar  in  order  to 
determine  how  much  radome  divorce  can  be  tolerated  and  still  maintain  a  desired  level  of 
tracking  accuracy. 
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L  DESCRIPTION 

The  FORTRAN  program  RADOME2  is  a  transmission-only  geometrical  optics  model 
of  an  antenna/radome  combination.  Its  main  purpose  is  to  predict  the  cross-polarization  that  is 
generated  when  the  radome  has  a  significant  level  of  divorce.  Co-pol  and  cross-pol  patterns  are 
computed  in  each  of  three  planes:  azimuth,  elevation  and  diagonal.  Monopulse  difference 
patterns  (elevation  or  azimuth)  may  also  be  calculated.  The  radome  shape  may  be  conical  or 
tangent  ogive,  either  one  with  user-specified  base  radius  and  height.  Alternatively,  the  radome 
may  be  removed  entirely  to  observe  the  antenna  properties  alone.  The  anienna  is  a  flat  plate 
array  with  a  circular  aperture,  and  may  be  scanned  mechanically  or  electronically.  It 
incorporates  a  low-sidelobe  Taylor  amplitude  distribution. 

The  code  is  not  intended  to  provide  accurate  estimates  of  sidelobe  levels.  When  the 
radome  has  significant  transmission  loss,  internal  reflections  will  be  significant,  and  may 
dominate  the  sidelobe  levels.  Those  internal  reflections  are  not  accounted  for  in  RADOME2. 
See  section  3  of  this  report  for  a  detailed  description  of  the  theory  and  derivation  for  the 
geometrical  optics  analysis.  This  appendix  is  intended  to  provide  instructions  for  using  the  code, 
as  well  as  enough  details  of  its  structure  to  allow  for  future  modifications. 

1.1.  Radome  Description 

The  user  may  choose  between  two  radome  geometries:  a  right  circular  cone;  and  a 
tangent  ogive.  In  both  cases,  only  two  independent  variables  need  to  be  specified:  the  base 
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radius;  and  the  "height,"  the  length  from  base  to  tip.  One  may  also  choose  no  radome,  in  which 
case  the  antenna  properties  alone  are  computed. 

1.2.  Antenna  Desci  ption 

The  antenna  is  a  circular  aperture,  whose  radius  is  less  than  that  of  the  radome.  Its 
center  is  always  located  at  (x,y,z)  =  (0,0,0),  which  is  also  the  center  of  the  radome  base.  The 
antenna  elements  are  spaced  in  an  equilateral  triangular  lattice.  The  maximum  number  of 
elements  is  500.  If  the  antenna  is  to  be  mechanically  scanned,  the  element  spacing  may  be 
chosen  as  anything  less  than  1  wavelength.  If  it  is  to  be  electronically  scanned,  the  elements 
should  be  closer  together  to  prevent  grating  lobes.  For  example,  a  spacing  of  about  .6X  is 
required  to  prevent  grating  lobes  at  scan  angles  of  60°. 

The  amplitude  weighting  is  a  Taylor  distribution,  which  gives  low  sidelobes.  The  user 
must  specify  two  parameters:  SLL,  the  dB  level  of  the  highest  sidelobe  relative  to  the  main 
beam;  and  NBAR,  an  integer,  which  is  the  number  of  sidelobes  at  the  peak  level.  The  pattern 
will  have  NBAR  sidelobes  on  each  side  of  the  main  beam  at  -SLL  dB,  and  the  remaining 
sidelobes  will  be  progressively  lower. 

The  user  must  specify  whether  the  antenna  is  to  be  mechanically  scanned,  or 
electronically  scanned.  If  electronic  scanning  is  chosen,  the  main  beam  will  be  somewhat 
broader  in  the  direction  of  scan.  In  either  case,  the  antenna  patterns  are  calculated  relative  to 
the  antenna’s  surface  normal.  The  antenna  is  y  polarized,  and  the  "azimuth"  and  "elevation" 
planes  are  taken  as  the  x-z,  and  y-z  planes,  respectively. 

Monopulse  difference  patterns  may  also  be  computed.  The  difference  pattern  is  generated 
by  phase  shifting  one  side  of  the  aperture  by  180°  relative  to  the  other.  The  usei  may  choose 
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either  elevation  difference  or  azimuth  difference.  Note  that  the  program  must  be  executed  three 
times  to  generate  all  three  monopulse  components:  E  (sum);  Aaz;  and  AEL. 

2.  INPUTS  AND  OUTPUTS 

RADOME2  reads  two  data  files  and  generates  four  others.  This  section  provides  a 
detailed  description  of  each  file.  Note  that  each  of  the  output  files  is  overwritten  if  it  already 
exists. 

The  first  input  file,  RADOME2.INP,  contains  information  on  the  geometry  and  other 
problem  parameters  such  as  scan  angles,  user  options,  etc.  It  is  an  optional  file:  the  user  may 
choose  instead  to  enter  all  the  same  data  from  the  keyboard.  The  second  input  file  contains  the 
transmissivity  data,  and  the  filename  is  entered  either  from  RADOME2.1NP  or  from  the 
keyboard. 

The  output  file  RADOME2.LOG  is  a  lengthly  "log"  file  containing  calculations  from 
each  stage  of  the  radome  transmission  calculation.  The  user  specifies  whether  or  not  the  file  is 
to  be  written.  Since  it  is  intended  primarily  for  validation  and  troubleshooting,  it  is  normally 
not  created.  RADOME2.PAT  is  the  output  file  for  antenna  radiation  patterns.  Its  seven 
columns  are  the  angle  and  co-  and  cross-pol  components  in  each  of  three  planes:  azimuth, 
elevation  and  diagonal.  RADOME2.APR  is  the  "equivalent"  aperture  distribution,  i.e.  the 
original  Taylor  distribution  modified  by  the  radome  transmission  properties.  It  contains  the 
element  locations  (original,  unscanned  x  an.'  y  coordinates)  and  the  co-  and  cross-pol  element 
excitations  as  complex  numbers.  RADOME2.IPT  contains  the  coordinates  of  where  the  ray 
from  each  antenna  element  intersects  the  radome  surface.  Finally,  RADOME2.SNM  contains 
the  three  components  of  the  surface  normal  at  the  intersect  point. 
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2.1.  Geometry  &  Configuration  File  (Input) 

The  first  input  the  program  requires  is  whether  the  configuration  data  is  to  be  entered 
from  the  keyboard  or  from  a  data  file.  If  the  user  responds  with  1  for  data  file,  the  program 
will  retrieve  all  remaining  inputs  from  the  file  RADOME2.INP.  An  example  of  this  file  is  given 
below.  Note  that  if  keyboard  input  is  chosen,  the  sequence  of  inputs  is  identical  to  those  in  the 
RADOME2.INP.  File  input  is  recommended,  since  most  of  the  parameters  remain  unchanged 
from  run  to  run. 


In  the  input  file,  there  may  be  any  number  of  comment  lines  preceding  the  data,  and  the 
last  comment  line  must  have  the  character  0  in  column  1.  If  there  are  no  comments,  that  will 
be  the  first  line.  The  data  is  unformatted,  but  the  values  must  appear  in  the  order  shown  in  the 
example  below. 


Input  file  for  program  RADOHE2 
Linel t  1  denotes  ogive  radome,  0  denotes 
cone  and  99  is  no  radome 
Line2i  7 . S  and  30.0  are  the  radome  base  radius 
and  height,  respectively,  in  wavelengths 
Line3i  5  -  antenna  aperture  radius, 

.5  »  interelement  spacing  (wavelengths) 

40  *  sidelobe  level 
4  -  NEAR 

Line4:  30.0,0.0  •  thetaO,  phiO 
Lines i  0PT2=0  for  mechanical  beam  eteering,  1  for  electronic 

OPT3*=0  for  sum  patterns,  1  for  az.  diff.,  2  for  el.  dlff. 
OPT4-1  enable  log  file,  anything  else  disables  logging 
Line6s  Data  file  for  transmissivity  values 
/  (last  comment  line  must  begin  with  this  character) 

1 

7.5,  30 

5.0,  0.5,  40,  4  Input  1b  unformatted 

30.0,  0.0 

0,  0,  1 

divl.dat 


Figure  Al.  Example  Input  File  RADOMH2.1NP 


The  scan  angles  are  are  spherical  coordinate  angles,  with  0  measured  from  the  z 


The  file  may 
contain  aa  many 
comment  lines 
aa  desired,  but 
they  must  all 
precede  the  data. 
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axis  and  4>  measured  from  the  x  axis  in  the  x-y  plane.  Thus  <f>0  is  the  scan  plane,  where  0°  is 
azimuth  and  90°  is  elevation. 


2.2.  Transmissivity  Data  File  (Input) 

An  example  transmissivity  data  file  is  shown  below.  Like  the  file  RADOME2.INP,  the 


data  is  preceded  by  a  line  with  ti  in  the  first  column,  and  all  lines  preceding  it  are  ignored. 


File 

DIV1 

DAT 

Sample 

file 

with  1. 

0  dB  and  10  deg 

divorce 

at  60 

degrees  incidence 

The  tile  may  contain  as  many 

transmissivity  vs. 

angle 

comment  lines  as  desired,  but 

they  must  all  precede  the  data. 

Para 

Perp 

Angle 

* 

Mag 

Phase 

Mag 

Phase 

(last  comment  line 

must  begin  with  §) 

19 

This  la  the  number  of 

values  in  the  file 

0.1 

0.0 

-11.4 

0.0 

-11.4 

S.O 

0.0 

-11.5 

0.0 

-11.2 

Input  data  is 

10. 0 

0.0 

-11.6 

0.0 

-11.1 

unformatted 

15.0 

0.0 

-11.7 

0.0 

-10.9 

20.0 

0.0 

-11.9 

0.0 

-10.6 

Magnitude  values 

25. C 

0.0 

-12.1 

0.0 

-10.1 

are  in  dB 

30.0 

0.0 

-12.5 

-0.0 

-9.7 

35.0 

0.0 

-12.9 

-0.0 

-9.1 

Phase  values 

40.0 

-0.1 

-13.4 

-0.0 

-8.4 

are  in  degrees 

45.0 

-0.2 

-13.9 

-0.0 

-7.8 

50.0 

-0.4 

-14.5 

-0.0 

-7.1 

55.0 

-0.7 

-15.2 

-0.0 

-6.6 

60.0 

-1.1 

-15.7 

-0.0 

-5.8 

65.0 

-1.6 

-16.3 

-0.0 

-4.6 

70.0 

-2.8 

-16.7 

0.0 

-3.7 

a 

75.0 

-4.2 

-16.7 

0.0 

-2.8 

80.0 

-7.5 

-16.5 

0.0 

-1.9 

85.0 

-15.0 

-15.8 

0.0 

-1.0 

90.0 

-45.0 

-14.7 

0.0 

-0.0 

Figure  A2.  Example  Transmissivity  Input  File 


The  first  data  line  is  the  number  of  angles  that  are  tabulated  in  the  file.  For  each  angle, 
the  program  reads  the  angle  in  degrees,  and  the  transmissivity  values  as  magnitude  in  dB  and 
phase  in  degrees. 
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The  program  will  use  this  input  data  to  create  lookup  tables  for  parallel  and  perpendicular 
transmissivities  at  1°  increments  from  0°  to  90°.  Hence  the  input  data  may  be  coarsely  sampled 
as  long  as  transmissivity  is  a  fairly  smooth  function  of  angle.  The  program  uses  Lagrange 
interpolation  to  create  the  lookup  table  from  the  input  data.  It  will  in  fact  extrapolate  if  the  input 
data  does  not  cover  the  entire  range  from  0°  to  90°.  Note  that  the  radomc  material  is  assumed 
to  be  anisotropic,  implying  that  its  transmissivity  is  a  function  only  of  the  angle  from  the 
normal. 


2.3.  Aperture  Distribution  File  (Output) 


El ement 

Coordinates 

Element 

Excitation 

Number 

in  wavelengths* 

Co-polarized 

Cross- 

polarized 

* 

X 

y 

real 

imag 

real 

imag 

1 

-.2500 

-4.9796 

.1164 

-.0146 

.1164 

-.0088 

2 

.2500 

-4.9796 

.1134 

-.0134 

.1134 

-.0083 

3 

-2.0000 

-4.5466 

.1223 

-.0173 

.1223 

-.0070 

4 

-1.5000 

-4.5466 

.1254 

-.0168 

.1254 

-.0076 

S 

-1.0000 

-4.5466 

.1312 

-.0169 

.1312 

-.0085 

6 

-.5000 

-4.5466 

.1354 

-.0162 

.1354 

-.0094 

7 

.0000 

-4.5466 

.1351 

-.0151 

.1351 

-.0092 

8 

.5000 

-4.5466 

.1303 

-.0134 

.1303 

-.0086 

9 

1.0000 

-4.5466 

.1223 

-.0114 

.1223 

-.0078 

10 

1.5000 

-4.5466 

.1105 

-.0083 

.1105 

-.0057 

• 

• 

one  line 

per  array 

element 

356 

-.7500 

4.5466 

.1336 

-.0163 

.1336 

.0091 

357 

-.2500 

4.5466 

.1362 

-.0160 

.1362 

.0095 

358 

.2500 

4.5466 

.  1335 

-.0146 

.1335 

.0090 

359 

.7500 

4.5466 

.  1270 

-.0127 

.1270 

.0083 

360 

1.2500 

4.5466 

.  1166 

-.0099 

.1166 

.0068 

361 

1.7500 

4.5466 

.1047 

-.0066 

.1047 

.0044 

362 

.0000 

4.9796 

.1144 

-.0138 

.1144 

.0084 

Figure  A3.  Example  Output  File  RADOME2.APR 


The  figure  above  shows  the  content  and  format  of  the  aperture  distribution  output  file. 
The  first  column  is  the  element  number  identifier.  The  next  two  columns  are  the  x,y  coordinates 
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of  the  unscanned  antenna  elements.  The  remaining  columns  are  the  effective  element  excitations 
accounting  for  the  transmission  properties  of  the  radome.  The  values  in  this  figure  are  those 
generated  using  the  input  files  RADOME2.INP  and  DIV1.DAT  shown  earlier. 


2.4.  Far  Field  File  (Ouput) 


Observation 

.Relative 

Rover  in 

dB 

Angle 

Azimuth 

Diagonal 

Elevation  j 

(deg) 

Co-pol 

X-pol 

Co-pol 

X-pol 

Co-pol 

X-pol 

-90 

-100.000 

-100.000 

-100.000 

-100.000 

-100.000 

-100.000 

-89 

-98.508 

-100.000 

-89.059 

-100.000 

-92.210 

-100.000 

-88 

-92.239 

-100.000 

-83.002 

-100.000 

-86. 185 

-100.000 

-87 

-88.326 

-98.287 

-79.422 

-98.014 

-82.660 

-100.000 

-86 

-85.319 

-95.843 

-76.847 

-95.623 

-80.161 

-100.000 

-85 

-82.790 

-93.976 

-74.821 

-93.819 

-78.233 

-100.000 

-10 

-40.983 

-90.525 

-38.877 

-42  o 453 

-37.424 

-46.419 

-9 

-31.418 

-85.915 

-30.407 

-39.914 

-29.411 

-41.904 

-8 

-25.167 

-82.909 

-24.557 

-37.992 

-23.838 

-38.782 

-7 

-20.537 

-80.283 

-20.138 

-36. 688 

-19.600 

-36.649 

-6 

-16.953 

-77.385 

-16.685 

-36.006 

-16.277 

-35.339 

-5 

-14.148 

-74.536 

-13.966 

-35.976 

-13.659 

-34.800 

-4 

-11.972 

-72.166 

-11.850 

-36.681 

-11.625 

-35.068 

-3 

-10.337 

-70.443 

-10.258 

-38.304 

-10.102 

-36.302 

-2 

-9.185 

-69.390 

-9.140 

-41.288 

-9.042 

-38.929 

-1 

-8.483 

-68.985 

-8.463 

-47.062 

-8.417 

-44.381 

0 

-8.210 

-69.195 

-8.210 

-69.195 

-8.210 

-69.195 

1 

-8.359 

-69.969 

-8.375 

-47.930 

-8.417 

-44.515 

2 

-8.933 

-71.202 

-8.963 

-42.734 

-9.042 

-38.985 

• 

85 

-81.854 

-93.488 

-74.674 

-92.974 

-78.370 

-100.000 

86 

-84.256 

-95.402 

-76.676 

-94.749 

-80.297 

-100.000 

87 

-87.150 

-97.884 

-79.231 

-97.130 

-82.796 

-100.000 

88 

-90.975 

-100.000 

-82.798 

-100.000 

-86.321 

-100.000 

89 

-97.186 

-100.000 

-88.847 

-100.000 

-92.346 

-100.000 

90 

-100.000 

-100.000 

-100.000 

-100.000 

-100.000 

-100.000 

Figure  A4.  Example  Pattern  Output  File  RADOME2.PAT 


The  listing  shown  above  is  the  pattern  file  generated  using  the  input  files  RADOME2.1NP 
and  DIV1.DAT  shown  earlier.  The  number  of  pattern  values  and  the  angle  increment  are  fixed 
by  a  PARAMETER  statement  in  the  main  program.  The  observation  angle  is  always  relative 
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to  the  antenna’s  surface  normal,  so  the  main  beam  appears  at  0°,  independent  of  the  scan  angle. 
The  values  in  RADOME2.PAT  may  represent  E  (sum)  Aaz  or  AEL,  depending  on  the  choice 
of  OPT(3)  in  the  input  file.  The  figure  above  shows  sum  patterns.  Values  are  in  dB  relative 
to  a  uniformly  illuminated  aperture.  Thus,  "taper  loss"  (the  reduction  in  aperture  efficiency  due 
to  a  low-sidelobe  distribution)  and  radome  transmission  loss  are  both  accounted  for,  and  that  is 
the  reason  that  the  main  beam  peak  is  lower  than  0  dB.  Note  that  the  patterns  are  "clipped"  so 
that  values  below  -100  dB  are  set  to  -100  dB. 

2.5.  Intercept  Points  and  Surface  Normals  (Outputs) 

The  two  figures  on  the  following  page  are  example  outputs  files  RADOME2.SNM  and 
RADOME2.IPT,  which  are  the  surface  normals  and  intercept  points,  respectively,  of  rays 
originating  at  each  antenna  element  and  passing  through  the  radome  to  an  equivalent  aperture 
outside  the  radome.  The  surface  normals  are  unit  vectors,  directed  inward.  In  the  .SNM  file, 
the  second  and  third  columns  are  the  x,y  coordinates  of  the  antenna  elements  before  scanning. 
In  the  .IPT  file,  the  first  column  is  the  antenna  element  index,  and  the  remaining  three  are  the 
intercept  coordinates  in  units  of  wavelenghts. 

2.6.  Log  File  (Output,  Optional) 

An  example  of  the  log  file  RADOME2.LOG  is  shown  on  the  following  page.  It  contains 
detailed  information  about  the  radome  transmission  calculations  made  for  each  element.  The  last 
option  in  the  input  file  is  set  to  1  if  this  file  is  to  be  created  (if  it  does  not  already  exist).  The 
first  block  of  information  written  is  a  verification  of  the  configuration  and  geometry  parameters. 
After  that,  a  separate  block  is  written  for  each  element.  If  the  meaning  of  this  output  data  is  not 
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Element 

Element 

Coordinatoa 

Components  of 

Unit 

Number 

in  Wavelengths* 

Surface 

Normal 

Vector 

# 

X 

y 

X 

y 

z 

1 

-.2500 

-4.9796 

-.4908 

.5124 

.7046 

2 

.2500 

-4.9796 

-.5242 

.5313 

.6655 

3 

-2.0000 

-4.5466 

-.3963 

.4002 

.8263 

4 

-1.5000 

-4.5466 

-.4273 

.4151 

.8032 

5 

-1.0000 

-4.5466 

-.4593 

.4308 

.7768 

6 

* 

-.5000 

-4.5466 

-.4923 

.4471 

.7468 

• 

359 

.7500 

4.5466 

-.5776 

-.4905 

.6525 

360 

1.2500 

4.5466 

-.6123 

-.5085 

.6055 

361 

1.7500 

4.5466 

-.6464 

-.5264 

.5523 

362 

.0000 

4.9796 

-.5074 

-.5218 

.6857 

Figure  A5.  Example  Output  File  RADOME2.SNM  (Surface  Normals) 


Aperture 

Ray-Radome 

Intercept 

Coordinates 

Element 

(in 

w avelengtha ) 

Number 

X 

y 

z 

1 

4.7693 

-4.9796 

8.7607 

2 

4.9125 

-4.9796 

8.0087 

3 

4.5024 

-4.5466 

11.7984 

4 

4.6802 

-4.5466 

11.1064 

5 

4.8480 

-4.5466 

10.3970 

6 

• 

5.0057 

-4.5466 

9.6701 

• 

359 

5.3543 

4.5466 

7.7739 

360 

5.4747 

4.5466 

6.9824 

361 

5.5836 

4.5466 

6.1711 

362 

4.8423 

4.9796 

8.3872 

Figure  A6.  Example  Output  File  RAD0ME2.IPT  (Intercept  Points) 


self-evident,  see  Section  3  of  this  report  for  a  detailed  explanation. 


The  figure  below  shows  the  order  in  which  the  program  executes.  The  boldface  titles 
next  to  flowchart  blocks  indicate  the  important  subroutines  that  assist  those  tasks.  The  most 
analytically  detailed  procedures  are  the  top  three  blocks  in  the  right  column.  Two  of  those 
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RAD0ME2  -  LOG  FILE 
ANTENNA  PARAMETERS: 

APERTURE  RADIUS  -  5.00 

SIDELOBE  LEVEL  *  -  40.  dB  NBAR  *  4 

SCAN  ANGLES:  THETA  *  30.00  dog.  PHI  -  .00  deg. 

RADOME  PARAMETERS: 

BASE  RADIUS  -7.50  HEIGHT  -30.00 
OPTIONS:  #1  #2  #3  #4 

10  0  1 


POINT  #  1  INTERSECT: 

SURFACE  NORMAL: 
INCIDENCE  ANGLE: 
E-CROSS-K: 
U-PERP: 
U-PARA: 

Incident  E  components 
Incident  H  components 
Co-pol  trans.  coeff. 

X  -pol  trans.  coeff. 


4.769  -4.980  8.761 

-.491  .512  .705 

68.601940 

-.866  .000  .500 

-.477  -.835  .275 

-.723  .550  .417 

(para,perp)t  .5504  -.8349 
(para,perp):  .8349  .5504 

(real,imag):  .7422  -.2206 

j real, imag) :  .9975  -.0702 


The  above  10  line3  are  repeated  for  each  antenna  element 


Figure  A7.  Example  Log  File  RADOME2.LOG 


calculations  are  bypassed  if  there  is  no  radome  present.  Further  details  of  the  actions  taken 
within  those  blocks,  as  well  as  within  the  two  subroutines  RCCONE  and  TANGIVE,  may  be 
found  in  Section  3  of  this  report,  as  well  as  in  the  internal  comments  in  the  source  listing  that 
follows. 
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PROGRAM 
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APPLICABLE 

SUBROUTINE 


RCCONE 

or 

TANGIVE 


PATTERN 


Figure  A8.  Flowchart  for  Program  RADOME2 
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4.  80URCB  CODE  LISTING 


RADOME 2  VI.  0 


Dan  McGrath 

Air  Force  Institute  of  Technology 

-  Department  of  Electrical  and  Computer  Engineering 

Wright-Pattereon  AFB,  OH  45433 

-  Origination  Date/  Last  modified!  25  Jul  91  /  15  Oct  91 
~  Language/  Operating  Syetemi  FORTRAN  (Microsoft)/  MS-DOS 


Purpose: 

This  program  simulates  the  pattern  performance  of  an  array 
antenna  with  a  non-ideal  radome.  The  antenna  may  be  either  phase 
steered  or  mechanically  scanned.  The  radome  geometry  may  be  either 
a  right-circular  cone  or  a  tangent  ogive.  The  output  consists  of 
the  co-  and  cross-polarized  azimuth  and  elevation  radiation  patterns. 
An  optional  output  is  the  co-  or  cross-polarized  aperture  amplitude 
distribution.  The  antenna  has  a  low-sidelobe  Taylor  distribution. 

The  radome  may  have  a  variable  transmissivity  vs.  incidence  angle, 
which  may  be  different  in  the  parallel  and  perpendicular  polari¬ 
zations. 

Input  Files: 

The  user  must  specify  an  Input  file  containing  the 
radome 's  parallel  and  perpendicular  plane  transmissivities 
(in  dB)  versus  incidence  angle  (in  deg.).  See  the  header 
for  routine  INPUT2  for  a  detailed  description  of  the  file 
format . 

The  problem  parameters  may  be  input  from  a  data  file 
or  manually  from  the  keyboard.  If  the  former  option  is 
chosen,  then  the  file  RADOME2 . INP  must  be  in  the  current 
directory.  See  subroutine  INPUT1. 

Output  Files: 

The  ouput  file  RAD0ME2.APR  will  contain  the  effec¬ 
tive  aperture  distributions:  the  antenna  element's 
coordinates  and  the  co-  or  cross-  polarized  amplitude 
and  phase  (in  degrees).  See  the  header  for  routine 
OUTPUT2  for  a  detailed  description  of  the  file  format. 

The  file  RADOME 2 . PAT  will  contain  all  four 
radiation  patterns:  co-  and  cross-pol  azimuth  and 
elevation  patterns  in  dBi. 

The  files  RADOME2 . IPT  and  RADOME2 . SNM  will  contain, 
respectively,  the  x,y,z  coordinates  and  the  surface  normal 
vector  components  at  the  points  where  rays  from  the 
antenna  intersect  the  radome  surface. 

All  four  of  the  above  output  files,  RADOME2.APR, 
RAD0ME2 . PAT ,  RADOME2 . IPT  and  RADOME2 . SNM  will  be  over¬ 
written  if  they  already  exist. 

Subroutines: 

CTAYLOR  -  Calculates  a  Taylor  amplitude  distribution  for 
a  circular-aperture  antenna 

BESSELJ  -  Function  subprogram  for  Bessel  functions  needed 
for  CTAYLOR 

PATTERN  -  Calculates  far  field  radiation  pattern  for  an 
array  antenna 

ANGLE  -  Finds  the  polar  angle  of  an  x,y  point 

INPUT1  -  Reads  problem  parameters 

INPUT2  -  Reads  radome  transmissivity  file 
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TANGIVE  -  Calculates  intersection  point  and  surface 
normal  vector  for  a  tangent  ogive 
RCCONE  -  Calculates  intersection  point  and  surface 
normal  vector  for  a  right  circular  cone 
PHASE  -  Finds  phase  angle  of  a  complex  number 
ATANH  -  Calculates  hyperbolic  arctangent  (for  CTAYLOR) 


PROGRAM  RAD0ME2 


Variable  Declarations 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


+ 


+ 


PARAMETER  (MXEL»500,MXRW«41,MXPF-361,MXNTA«90) 

MX  EL  -  Maximum  number  of  antenna  array  elements 
MXRW  »  Maximum  number  of  rows  in  the  array 
MXFF  »  Maximum  number  of  far  field  pattern  angles 
MXNTA=  Max  number  of  angles  in  transmissivity  tables 
PARAMETER  (PI=3. 14159265359, DTR-P1/180. ) 

DTP,  *  Degrees  to  radians  conversion 
PARAMETER  ( PT0=-90 DT-1 . 0 , NANG* 181 ) 

PTO  *»  Initial  theta  angle  for  pattern  calculations 
DT  ■  theta  increment 
NANG  *  number  of  pattern  points 
DIMENSION  X (MXEL) ,Y (MXEL) , Z (MXEL) , R (MXEL) , AMP (MXEL) , PHS (MXEL) 
X,Y,Z  =  Antenna  element  coordinates 
(dimensions  in  wavelengths) 

R  =  Radius  to  element  from  array  center  (wavelengths) 

AMP, PHS  «  Scratch  arrays  for  element  amplitude  and 
phase  for  pattern  calculation 
COMPLEX  A(MXEL) , AC(MXEL) , AX(MXEL) 

A  =■  Array  element  excitation 

AC, AX  *  Effective  co-pol  and  cross-pol  aperture 
amplitudes  accounting  for  radome  effects 
DIMENSION  XIP(MXEL) , YIP(MXEL) , ZIP(MXEL) , 

XSN(MXEL) , YSN(MXEL) , ZSN(MXEL) 

XIP,  YIP,  ZIP  »■  Coordinates  of  the  point  where  a  ray 

from  the  array  element  intercepts  the  radome 
XSN, YSN, ZSN  *  Surface  normal  vector  component*  at 
the  intercept  point 

DIMENSION  FFCA(MXFF) , FFCE (MXFF ), FFCD ( MXFF) , 

FFXA ( MXFF ), FFXE ( MXFF ), FFXD ( MXFF ) 

FFCA  «=  Co-polarized  azimuth  pattern  in  dBi 
FFCE  =  Co-polarized  elevetion  pattern 
FFCD  =  Co-polarized  diagonal  plane  pattern 
FFXA  «*  Cross-polarized  azimuth  pattern 
FFXE  ■  Cross-polarized  elevation  pattern 
FFXD  ■  Cross-polarized  diagonal  pattern 
COMPLEX  TCA ( 0 : MXNTA ) , TXA ( 0 : MXNTA ) 

TCA  ■  Co-pol  transmissivity  amplitude  (dB) 

TXA  »  Cross-pol  transmissivity  amplitude  (dB) 

COMPLEX  ETL, ETM 


ETL,ETM  are  temporary  variables  used  for 

transmissivity  calculations 
INTEGER  OPT (4 ) , NBAR 

OPT ( 1 )  *  0  for  conical  radome,  1  for  tangent  ogive 
or  99  for  no  radome 

OPT ( 2 )  *  0  for  mechanical  beam  steering,  1  for 
phase  steering 

OPT(3)  *  0  to  simulate  sum  channel  patterns 
*  1  for  azimuth  difference 
=  2  for  elevation  difference 

OPT ( 4 )  =  1  to  save  intermediate  calculations  in  file 
RADOME 2 . LOG 

NBAR  =  Number  of  fixed-level  sidelobes  in  antenna 

pattern  as  determined  by  Taylor  distribution 
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C 


CHARACTER* 20  FIL1 

PIL1  ■  Input  file  for  radome  transmissivity 
REAL  RAD , TO , P0 , BAS , HGT , SLL 

RAD  *  Antenna  radius  in  wavelengths 

TO, PO  -  Spherical  coordinate  angles  theta  and  phi 

for  scan  direction  {radome  axis  is  theta«0) 
in  radians 


BAS  -  Radome  base  radius  in  wavelengths 
HGT  m  Radome  height  in  wavelengths 

DD  -  Array  element  spacing  in  wavelengths 

SLL  =  Antenna  design  sidelobo  level  in  dB 
relative  to  main  beam  peak 


Input  and  Array  Element  Calculations 


Introduction  message 
WRITS ( * , 1 ) 

FORMAT (//,1X,30( -  PROGRAM  RADOME -2  - 

+  IX, 30 ('-'), //) 

Get  problem  parameters 

CALL  INPUT1 ( RAD, TO , P0 , FIL) , BAS , HGT, DD , SLL, NBAR, OPT) 

Get  radome  transmissivity  values  from  data  file 
and  interpolate  to  create  a  lookup  table 

IF(OPT(l) .NE.99)  CALL  INPUT2 {FILl, TCA, TXA) 

Calculate  the  element  locations  for  the  unscanned  antenna: 

These  are  an  equilateral-triangular  lattice  with 
spacing  DD  in  wavelengths  between  elements. 
DD«.5/sin(60)  is  adequate  to  prevent  grating  lobes 
for  a  phase-steered  array,  and  DD=l/ein(60)  is 
adequate  for  a  mechanically-steered  array. 

SIN60=SIN(PI/3. ) 

NEL=0 

DY=DD*SIN60 

XX=0 

NX=RAD/DD+5 
NY=NX/SIN60 
DO  10  I  =  -NY , NY+1 
XX=XX+ . 5 

IF (XX . GE . . 75 )  XX=XX-1. 

Yl* ( I- . 5 ) *DY 
DO  10  K=-NX , NX 
Xl=(K+XX) *DD 
RR3SQRT( XI* *2 +Y 1**2) 

IF(RR.GT.RAD)  GO  TO  10 
NEL=NEL+1 

IF(NEL.GT.MXEL)  THEN 

WRITE ( * , * )  '  >»  Error:  Too  many  array  elements' 

GO  TO  1111 
END  IF 
X ( NEL ) =X1 
Y ( NEL) =Y1 
Z ( NEL ) =0 . 0 
R(NEL) =RR 
CONTINUE 
WRITE(*, 12)  NEL 

FORMAT ( '  .  Number  of  array  elements®  ',14,'  . ') 
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C  Calculate  the  element  amplitude  weights 
C 

DO  20  I«1,NEL 
NCALL*I-1 
RNORM“R( I ) /RAD 

CALL  CTAY LOR { SLL, RNORM, NBAR, GP, NCALL) 

A  ( I  )  *=CMPLX  ( GP , 0 .  ) 

C  Reverse  amplitudes  as  appropriate  for 

C  azimuth  or  elevation  difference 

IF(OPT(3 ) . EQ. 1 )  THEN 

I F  (  X ( I )  .  LT .  0  .  )  A{  I )  *=CMPLX  ( -GP, 0 . ) 
IF(X(I) .EQ.O. )  A(I)-(0.,0.) 

ELSE  IF(OPT(3) .EQ.2)  THEN 
IF(Y(I) . LT. 0 . )  A { I ) =CMPLX ( -GP , 0 . ) 
IF(Y(I) .EQ.O. )  A(I)*=(0.,0.) 

END  IF 


C 


20 

C- 

c 

c- 

c 

c 

c 

c 

c 


41 


50 

c 

65 


91 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 


initialize  co  &  cross  pol  distributions 
AC( I ) =A( I ) 

AX (I)=(0. 0,0.0) 

CONTINUE 


Antenna  Scanning 


Calculate  the  array  element  phase  shifts  for  electronic 
scanning  to  the  spherical  angle  theta=T0,  phi=P0. 

Get  sin,  cos  of  scan  angles 
ST0=SIN (TO ) 

CTO=COS (TO ) 

SP0=SIN(P0) 

CP0=COS(P0) 

IF  (OPT(2) .EQ.O. OR. TO. EQ.O. )  GO  TO  65 
WK=2 . *PI 
WKS=WK*ST0 
DO  50  I«1,NEL 

PHAS=WKS* (X  ( I ) *CP0+Y ( I ) *SP0 ) 

AAMP=REAL  ( A ( I ) ) 

A ( I ) =CM?LX ( AAMP*COS ( PHAS ) , AAMP*SIN ( PHAS ) ) 

CONTINUE 

Begin  recording  to  log  file 
IF(OPT(4).EQ.l)  THEN 

OPEN (9, FILE= ' RADOHE2 . LOG ' ) 

WRITE ( 9 , 9 1 )  RAD , SLL ,  NBAR , T0/DTR , PO/DTR , BAS , HGT,  OPT 
FORMAT ( / ,  1 OX ,  '  RADOME2  -  LOG  FILE  ANTENNA 

+  •PARAMETERS:',/,'  APERTURE  RADIUS  - ' , F5 . 2 , / , ' SIDELOBE 
+  'LEVEL  -  - ' , F4 . 0, '  dB  NBAR  •',12,/,'  SCAN  ANGLES:  *, 

+  'THETA  =' ,F7.2, •  deg.  PHI  «',F7.2,’  deg.’,/,’  RADOME  ', 

+  'PARAMETERS: '  BASE  RADIUS  *’,F5.2,'  HEIGHT  «*,F5.2,/, 
+  '  OPTIONS:  #1  /2  #3  #4 ' , / , 12X, 12 , 3 ( 3X, 12 ) , / ) 

END  IF 

IF (OPT ( 1 ) . EQ . 99 )  GOTO  115 

WRITE ( * , * )  *  .  Beginning  Aperture  Calculation . ' 

WRITE (*,*) 


Projected  Aperture  Calculation 
If  there  is  no  radome,  thiB  block  is  skipped 


Calculate  the  effective  co-  and  cross-polarized  aperture 
amplitude  and  phase  distributions  as  follows: 

(1)  Find  intersection  of  ray  from  element  I  with  radome 

(2)  Find  radome  surface  normal  at  intersection  point 

(3)  Interpolate  from  radome  transmissivity  table  to 
obtain  para.  6  perp.  pol.  transm.  coefficients 

(4)  Apply  coefficients  to  decomposed  incident  wave 
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polarization,  than  recombine 


DO  100  1-1, NEL 

WRITE ( * , ' ( 1H+, 5X, A, 14 ) ' )  *  Element  number:  ',1 
In  the  case  of  mechanical  antenna  scanning, 
calculate  the  coordinates  of  the  gimballed 
antenna  array  element 
IF(OPT(2) . EQ.O. AND.T0.NE.0. )  THEN 

skip  this  calculation  if  the  scan  angle  TO-O 
XR*X ( I ) * (CTO*CPO**2+SPO*  *2 ) +Y ( I ) *CP0*SP0* ( 1-CT0) 

YR«X ( I ) *CP0*SP0* ( 1-CT0) +Y ( I )  * (CT0*SP0**2+CP0**2 ) 

ZR— X ( I ) *ST0*CP0-Y ( I ) *ST0*SP0 
ELSE 
XR“X( I ) 

YR»Y ( I ) 

2R«Z(I) 

END  IP 

Pind  coordinates  of  aperture  points 
XA“XR+HGT*ST0*CP0 
YA=YR+HG'  ' ST0*SP0 
ZA=ZR+1:  1  TTO 

Find  intersect  points  and  surface  normals 
IF (OPT ( 1 ) . EQ. 1 )  THEN 

CALL  TANGIVE (XR, YR, ZR,  XA, YA, ZA, BAS , HGT, 

XiP  ( I )  ,  Y'lP  ( I )  ,  ZIP(I)  ,  XSN  (  I )  ,YSN(I)  ,  ZSN  ( I )  ) 

ELSE 

CALL  RCCONE  ( XR , YR , ZR , XA , YA , 2 A , BAS ,  HGT , 

XIP(I) ,YIP(I),2IP(I) ,XSN(I),YSN(I),ZSN(I) ) 

END  IP 

I F ( OPT ( 4 ) . EQ . 1 )  THEN 

WRITE (9, 64)  I,XIP(I) ,YIP(I),ZIP(I),XSN(I) , YSN ( I ) , ZSN ( I ) 
FORMAT ( 1  POINT#', 13, •  INTERSECT: ’, 3P7 . 3 ,/ , 

'  SURFACE  NORMAL: ’ ,3F7. 3) 

END  IP 

Calculate  incident  polarization  ("e")  vector 
Note  that  the  antenna  is  assumed  to  be  y-pol- 
arized  before  scanning.  It  remains  so  if  the 
steering  is  electronic. 

I F ( OPT ( 2 ) . EQ . 0 )  THEN 
XP=CPO*SPO* (1-CT0) 

YP=CTO*SPO**2+CPO**2 

Z?=-ST0*SP0 

ELSE 

XP=0. 

YP=1 . 0 
ZP=0. 

END  IF 

Calculate  incidence  angle:  the  normalized 
dot  product  of  the  ray  direction  with  the 
surface  normal  is  the  cosine  of  the 
incidence  angle 

XRD, YRD, ZRD  are  the  "k"  vector 
RDMG=SQRT( ( XA-XR ) * *2+ ( Y A-YR ) *  * 2 + ( ZA- ZR ) * *2 ) 

XRD= (XA-XR) /RDMG 
YRD=(YA-YR)/RDMG 
ZRD= ( ZA-ZR) /RDMG 

AINC=ACOS (ABS(XSN(I)*XRD+YSN(I) * YRD+ZSN (I)*ZRD))*180./PI 
I INC=AINC 

IF ( I INC . GT. 90 )  I INC=180-IINC 

IF(OPT(4) .EQ. 1)  WRITE ( 9 , * )  '  INCIDENCE  ANGLE AINC 

Find  cross  product  e  x  k 
XEK= ( ZP* YRD-YP*  ZRD ) 

YEK=(XP*ZRD-ZP*XPD) 

ZEK-(YP*XRD-XP*YRD) 

Find  vector  normal  to  plane  of  incidence 


c 

c 


c 


c 

c 

c 

c 


66 


C 

C 


C 


C 

C 

C 

C 

C 

C 


71 


100 
c — 
c 

c— 

115 

C 

C 

C 

C 

C 

C 

C 

c 


(The  notation  convention  usee  laat  letter  L  or  M 
to  denote  perpendicular  and  parallel,  respectively) 
XUL- ( ZSN ( I ) *YRD-YSN( I ) *ZRD ) 

YUL«(XSN(I)*ZRD-ZSN(I)*XRD) 

ZUL- ( YSN ( I ) ‘XRD-XSN ( I ) * YRD ) 

and  normalize  it  (this  is  "u-perp") 
ULM=SQRT(XUL**2+YUL**2+ZUL**2) 

XUL-XUL/ULM 

YUL-YUL/ULM 

ZUL«=ZUL/ULM 

Find  vector  in  plane  of  incidence  ("u-para") 

(This  one  doesn't  need  normalization  because  it 
is  the  cross  product  of  unit  vectors,  and  so  it 
must  also  be  a  unit  vector) 

XUM* ( YUL*ZRD-ZUL* YRD ) 

YUM* ( ZUL*XRD-XUL*ZRD ) 

ZUM* ( XUL*YRD-YUL*XRD ) 

IF (OPT (4) .  EQ, 1)  THEN 

WRITE (9, 66)  XEK, YEK, ZEK,XUL,YUL, ZUL, XUM, YUM, ZUM 
FORMAT ( '  E-CROSS-K: ' , 3F7.3, /, 

+  '  U-PERP: ' , 3F7 . 3, /, 

+  '  U-PARA: ' ,3F7. 3) 

END  IF 

Decompose  incident  wave  into  parallel  and 
perpendicular  components 
EIL=(XP*XUL+YP*YUL+ZP*ZUL) 

EIM= ( XP*XUM+YP* YUM+ZP*  ZUM) 

HIL*  ( XEK*XUL+YEK*YUL+ZEK*ZUL) 

HIM* ( XEK*XUM+YEK*YUM+ZEK*ZUM) 

Calculate  transmitted  components 
ETL“TXA( I INC) *EIL 
ETM=TCA( I INC) *EIM 

AC ( I )  represents  the  dot  product  of  the  transmitted 
E  field  vector  with  the  unit  vector  in  the  incident 
E  field  direction 
AC( I ) =A( I ) * (ETL*EIL+ETM*EIM) 

AX ( I )  is  the  dot  product  of  the  transmitted  field 
E  field  vector  with  the  unit  vector  in  the  incident 
H  field  direction 
AX ( I ) =A( I ) * ( ET-*HIL+ETM*HIM) 

IF(OPT(4) . EQ. 1)  THEN 

WRITE(9, 71)  EIM, EIL, HIM, HIL, TCA ( IINC) , TXA( I INC) 

FORMAT ( '  Incident  E  components  (para,perp) : ' ,2F7.4, /, 

+  '  Incident  H  components  (para, perp) i ' , 2F7 . 4, / , 

+  '  Co-pol  trans.  coeff.  { real , imag ) : ' , 2F7 .4, / , 

+  '  X  -pol  trans.  coeff.  ( real, imag ) : ' , 2F7 . 4, / ) 

END  IF 
CONTINUE 


Radiation  Patterns 


WRITE ( * , * )  '  .  Calculating  Patterns . ’ 

The  pattern  cuts  must  be  made  relative  to  the  main 
beam  position.  In  the  phase  scanning  case,  it  is 
necessary  to  rot  ite  the  antenna  by  TO  in  the  (P0+180) 
direction  in  order  to  bring  the  main  beam  back  to  the 
theta, phi=0,0  direction: 

I F ( TO . EQ . 0 . . OR . OPT ( 2 ) . EQ . 0 )  GO  TO  119 

Skip  this  if  the  scan  direction  is  zero 
or  if  the  antenna  is  mechanically  scanned 
DO  117  1=1, NEL 

XR=X ( I ) * ( CT0*CP0*  *2  +  SPO**2)+Y(I) *CPO*SPO* ( 1-CT0) 

YR=X ( I ) *CP0*SP0* ( 1-CT0 ) +Y ( T ) * ( CTO * SFO* *2 +CP0* *2 ) 

ZR=X ( I ) *ST0*CP0+Y ; I ) *STO*SPO 
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X(I)-XR 
y(I)=yr 
Z(I)=ZR 
117  CONTINUE 
C  Co-pol  patterns 

119  DO  120  I-1,NEI. 

AMP  ( I )  *CABS  ( AC  ( I )  ) 

PHS ( I ) * PHASE (AC{ I )) *  180. /PI 

120  CONTINUE 

CALL  PATTERN ( X ,  Y ,  Z , AMP ,  PHS ,  NEL , PTO , DT , NANO ,  0.0, FFCA, 1.0,1) 
CALL  PATTERN (X,Y,Z, AMP, PHS, NEL, PTO, DT, NANG, 45. 0,PPCD, 1.0,1) 
CALL  PATTERN ( X ,  Y ,  Z , AMP , PHS , NEL , PTO , DT , NANO ,90.0, PFCE ,1.0,1) 

C  Crose-pol  patterns 

IF(OPT(l) .EQ.99)  GO  TO  132 
AMPMX«0. 

DO  130  1=1, NEL 

AMP( I ) “CABS ( AX( I ) ) 

I F ( AMP ( I ) . GT . AMPMX )  AMPMX-AMP(I) 

PHS ( I ) "PHASE (AX ( I))*1S0./?I 
130  CONTINUE 

IF (AMPMX. EQ.O)  THEN 

WRITE ( * , * )  1  lit  The  cross-pol  amplitude  is  zero;  cross-  ’ 
WRITE ( * ,  * )  '  pol  patterns  will  not  be  calculated.’ 

END  IF 

C  If  the  croaa-pol  amplitude  is  zero  or  if  there  is 

C  no  radome,  then  there  is  no  cross  polarized  pattern 

132  IF(AMPMX. EQ.O. OR. OPT(l) .EQ.99)  THEN 

DO  135  1*1, NANG 
FFXA(  I )  «=0 . 0 
FFXD  ( I )  *=0 . 0 
FFXE(I)=0.0 
135  CONTINUE 

ELSE 

CALL  PATTERNS,  Y,  Z, AMP, PHS, NEL,  PT0,DT, NANG,  0.0,  FFXA,  1 .0,  1 ) 
CALL  PATTERN (X,Y,Z, AMP, PHS, NEL, PTO, DT, NANG, 45.0, FFXD, 1.0, 1) 
CALL  PATTERN ( X , Y , Z , AMP , PHS , NEL , PTO , DT , NANG ,90.0, FFXE ,1.0,1) 
END  IF 

C - 

C  Output 

C - 

C 

IF(OPT(4) .EQ. 1)  CLOSE ( 9 ) 

C  Intersection  points 

OPEN ( 11 ,  FILE5*  '  RAUOME2  .  IPT '  ) 

DO  210  1*1 , NEL 

WRITE (11, 211)  I,XIP(I),YIP(I),ZIP(I) 

210  CONTINUE 

211  FORMAT (2X, I4,3(2X,F10.4) ) 

CLOSE (11) 

C 

c  Surface  normal  components 

OPEN ( 12 , FILE* ' RAD0ME2 . SNM* ) 

DO  220  1=1, NEL 

WRITE (12, 213)  I,X(I) ,Y(I) ,XSN(I) ,YSN(I) ,ZSN(I) 

213  FORMAT ( 2X ,I4,5(2X,F10.4)) 

220  CONTINUE 

CLOSE (12) 

C 

C  Aperture  distributions 

OPEN(13,FILE=' RADOME 2 . APR  '  ) 

DO  230  1=1, NEL 

IF(OPT(l) .EQ.99)  AX(I)=(0. 0,0.0) 

WRITE (13, 231)  I.X(I) ,Y(I) ,CAB3(AX( 1 ) ) , PHASE ( AX ( I ) ) , 

+  CABS ( AC ( I ) ) , PHASE ( AC ( I ) ) 

230  CONTINUE 
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231  FORMAT ( 2X, X4,6(2X,F10.5) ) 

CLOSE (13) 

C 

C  Far  field  patterns 

OPEN ( 14, FILE- 'RADOME2. PAT* ) 

DO  240  1=1 , NANO 
11=1-91 

WRITE (14 ,241)  II , FFCA ( I ) , FFXA( I ) , PFCD ( I ) , 
+  FFXD ( I ) , FFCE ( I ) , FFXE ( I ) 

240  CONTINUE 

241  FORMAT (2X, I4,6(2X,F10.5)) 

CLOSE (14) 

1111  STOP 
END 

C - 

C  End  of  main  program  RADOME2 

C - 


(||a=se  =  ssERSBSSSCsaauBiibSSeeBSEssss  =  ssBascss3cssssss  =  =  Bcszs9s&scn86sasBB 

C  SUBROUTINE  INPUT1 

C  This  routine  handles  all  input  and  error  checking  of  inputs.  It 
C  offers  the  user  the  option  of  entering  all  parameters  from  the 
C  keyboard,  or  from  the  data  file  RADOME2 . INP . 

C  The  input  file  format  is  as  follows: 

C  Comment  line3:  the  file  may  contain  a  header  with  any  number  of 
C  comment  lines,  the  last  of  which  contains  a  t  in  column  1. 

C  Data  records  are  in  free  format,  containing  the  following: 

C  1:  OPT ( 1 )  (Radome  type  option) 

C  2:  RAD, DO, SLL, NBAR  (Antenna  radius,  element  spacing, 

C  sidelobe  level  and  Taylor  NBATt  parameter) 

C  3:  TO, PO  (Scon  angles  theta, phi) 

C  4:  OPT ( 2 ) , OPT ( 3 ) , OPT ( 4 )  (Options) 

C  5:  FIL1  (Radome  transmissivity  data  file  if  0PT1<>99) 

^cccauEmBtnxnwKiiitDrsasDtiaasaiiaBDaxvBkattHBnnaaaBDBKaiizcBnsanascniBbCEszBSB 

SUBROUTINE  INPUT1 (RAD, TO,  PO, FIL1 , BAS, HOT, DD, SLL, NBAR, OPT ) 

INTEGER  OPT ( 4 ) 

CHARACTER*20  FIL1 
CHARACTER* 1  CHI 
PI-3.141592  359 

WRITE (*,* (A\ )' )  '  Enter  0  for  keyboard  input,  1  for  data  file:’ 
READ ( * , * )  I  KEY 
IF(IKEY.EQ.l)  GO  TO  51 
WRITE( * , 20) 

20  FORMAT (*  Choose  the  radome  type:',/,*  0  for  cone,’ 

+  ’  1  for  ogive  or  99  for  no  radome') 

READ ( * , * )  OPT ( 1 ) 

I F ( OPT ( 1 ) . NE . 9  9 )  THEN 
WRITE (*,21) 

21  FORMAT (iX,'  Enter  the  following  radome  parameters:',/, 

+  •  BAS  -  base  radiue  in  wavelengths',/,'  HGT  =  height', 

+  *  in  wavelengths') 

WRITE ( * , ' ( A\ ) ' )  '  Enter  BAS, HGT’ 

READ( * , *)  BAS, HGT 
END  IF 
WRITE ( * , 22  ) 

22  FORMAT (IX,  '  Enter  the  following  antenna  parameters:',/, 

+  ’  RAD  =  aperture  radius  in  wavelengths',/, 

+  '  DD  =  interelement  spacing  in  wavelenghts ' , / , 

+  '  SLL  -  peak  sidelobe  level,  e.g.  40',/, 

+  '  NBAR  -  Taylor  parameter',/,*  TO,PO  =', 

+  '  theta  and  phi  scan  angles  (deg)') 

WRITE ( * , ' ( A\ ) ' )  '  Enter  RAD, DD , SLL, NBAR, TO, P0:  ’ 
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READ(*,*)  RAD,DD,SLL,NBAR,TO,PO 
TO-TO*PI/180. 

PO-PO*PI/180. 

WRITE (*,23) 

23  FORMAT ( IX, 1  Choose  the  following  options:’,/, 

+  '  OPT2  -  0  for  mechanical  scanning,  or  1  for  electronic',/, 

+  '  OPT3  ■  0  for  sum  patterns,  ■!  for  az.  diff.,  "2', 

+  '  for  el.  diff.*,/,'  OPT4  •  1  to  save  intermediate', 

+  '  calculations  in  a  log  file') 

WRITE ( * ,  ' ( A\ ) ' )  *  Enter  0PT2,  0PT3,  0PT4:  ' 

READ ( * , * )  OPT ( 2 ) , OPT ( 3 ) , OPT ( 4 ) 

IF (OPT ( 1 ) . NE . 99 )  THEN 

WRITE ( * , * )  '  What  is  the  name  of  the  data  file  containing' 
WRITER,  *  (A\)  • )  '  the  transmissivity  data?  • 

READ ( * , ' ( A) ' )  FIL1 
END  IP 
RETURN 


C -  This  block  for  file  input  - - - 

51  OPEN(4  ,  FILE** '  RADOMB2  .  INP  '  ,  STATUS-'OLD*  ) 

C  Bypass  comment  lines  (if  any)  in  input  file 

55  READ ( 4 , ' (Al) ' )  CHI 


IF (CHI . NE . ' # ' )  GO  TO  55 
READ (4,*)  OPT ( 1 ) 

IF (OPT ( 1 ) . NE . SO )  READ ( 4 , * )  BAS, HOT 
READ ( 4 , * )  RAD, DD , SLL, NBAR 
READ (4,*)  TO, PO 
TO=TO*PI/180. 

PO=PO*PI/180. 

READ ( 4 , * )  OPT ( 2 ) , OPT ( 3 ) , OPT ( 4 ) 

IF (OPT (1) .NE.99)  READ (4, ' (A) ' )  FIL1 
RETURN 
END 


O^Ec  =  =  =  =  c::  =  cccEzc3BCDSXieKX=Ec:[:  =  cc3  =  acci:cDBzcxc  =  caBcc  =  rE:zc  =  Bzot  ==cc  =  =  =  =  z 

C  SUBROUTINE  INPUT2 

C  This  subroutine  reads  a  data  file  containing  radome  surface 
C  transmissivity  information  and  uses  interpolation  to  create  a 
C  lookup  table. 

C  The  input  data  file  may  contain  any  number  of  comment  lines 
C  preceding  the  actual  data.  The  symbol  #  muBt  be  in  the  first 
C  column  of  the  last  comment  line  (which  may  or  may  not  contain 
C  any  comment).  The  first  data  record  is  the  number  of  angles  for 
C  which  data  is  available,  NANG.  The  next  NANG  lines  contain  the 
C  incidence  angle  in  degrees,  the  co-pol  transmissivity  in  dB, 

C  the  co-pol  transmission  phase  in  deg,  the  croos-pol  transmissivity 
C  in  dB,  and  the  cross-pol  transmission  phase  in  deg.  Zero  degrees 
C  incidence  is  the  surface  normal  direction.  It  is  assumed  that 
C  the  radome  material  is  isotropic,  that  is  it  has  the  same  value 
C  of  transmissivity  at  the  angle  theta  in  any  phi  plane. 

C  Three-point  Lagrange  interpolation  is  used  to  fill  in  the 
C  lookup  table  entries  between  0  and  90  degrees.  0  degrees  is 

C  normal  incidence.  Note  that  Lagrange  interpolation  can  also  be 

C  used  to  extrapolate.  It  the  data  file  does  not  contain  values 
C  all  the  way  out  to  90  degrees  from  normal,  the  table  values  will 
C  be  extrapolated  out  to  90  degrees.  Unfortunately,  this  extrapo- 
C  lation  may  not  be  accurate,  especially  if  the  last  angle  given  is 
C  much  less  than  90  degrees.  Note  also  that  the  phases  in  the  input 

C  data  file  must  be  continuous,  without  "wrapping"  at  +/-180. 

C  The  output  data  is  complex  values  ("gamma"). 

QBaB  =  seexeE  =  =  ssEBESBBKBa  =  KeaBBS!csisussEsaEsscs  =  sz£2S£:t::s:n  =  scs3scss  =  sa 


SUBROUTINE  INPUT2 ( FIL1 , TC, TX ) 

COMPLEX  TC ( 91 ) , TX ( 9 1 ) 

C  TCA,TXA  =  co-  and  cross-pol  transmission  coefficients 

C  (sometimes  denoted  gamma-para/perp) 
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REAL  OCA ( 91 ) , DCP (91 ) , DXA ( 91 ) , DXP (91 ) , DANG (91 ) 

DCA, DCP, DXA, DXP  are  scratch  arrays  for  reading  the 
values  from  the  data  file.  DANG  is  an  array  with 
the  angle  values. 

CHARACTER* 20  PIL1 
CHARACTER* 1  CHI 

Read  the  data  file 
OPEN ( 10, FILE*FILl, STATUS* ' OLD  1 ) 

Bypass  comment  lines 
READ(10, ’ (Al) ' )  CHI 
IF(CH1.NE. ' )  GO  TO  5 

Get  the  number  of  data  points 
READ ( 10, * )  NANG 

Read  the  data 
DO  10  1*1 , NANG 

READ (10,*)  DANG ( 1 ) , DCA( I ) , DCP ( I ) , DXA ( I ) , DXP { I ) 
Convert  amplitudes  from  dB  to  volts 
DCA(I)*10.**(DCA(I)/20. ) 

DXA ( I ) *10 . ** (DXA{ I ) /20 . ) 

0  CONTINUE 
CLOSE (10) 

Begin  interpolation 

N=3 

DO  20  1*0,90 
J=I  +  1 
X-I 

IF(X. GT. DANG (N ) .AND. N.LT. NANG)  N=N+1 
X0*DANG (N-2 ) 

X1*DANG(N-1) 

X2«DANG(N) 

SO* ( X-Xl ) * (X-X2 ) / ( X0-X1 ) / ( X0-X2 ) 
S1»(X-X0)*(X-X2)/(X1-X0)/(X1-X2) 

S2=  ( X-X0 )  *  ^-Xl )  /  ( X2-X0 )  /  (X2-X1 ) 

TCA=DCA(N-2)*S0+DCA(N-1)*S1+DCA(N)*S2 
TCP*DCP (N-2 ) *S0+DCP ( N— 1 ) *S1+DCP ( N) *S2 
TXA=DXA(N-2 ) *S0+DXA(N-1 ) *S1+DXA(N) *S2 
TXP=DXP (N-2 ) *S0+DXP(N-1) *S1+DXP(N) *S2 
Convert  phase  to  radians 
DTR=3. 141592  359/180. 

TCP=TCP*DTR 

TXP=TXP*DTR 

Save  as  complex 

TC( J ) *CMPLX (TCA*COS (TCP ) ,TCA*SIN(TCP) ) 

TX ( J ) =CMPLX (TXA*COS (TXP ) ,TXA*SIN(TXP) ) 

0  CONTINUE 
RETURN 
END 


t=====e====EC=5£tCS=E=3K3tBC! 


SUBROUTINE  TANGIVE 

This  subroutine  finds  the  point  of  intersection  of  a  ray 
with  a  tangent  ogive  radome  and  the  surface  normal  vector  at  the 
point  of  intersection. 

The  inputs  are  XO,YO,ZO,  the  coordinates  of  the  ray's  initial 
point,  XA,YA,ZA,  the  coordinates  of  the  ray's  destination  point 
on  the  effective  aperture  outside  the  radome,  and  the  radome  base 
(BB)  and  height  (H) . 

The  point  where  the  ray  intersects  the  radome  is  XI,YI,ZI. 

The  components  of  the  surface  normal  vector  are  XSN,YSN,ZSN. 

All  dimensions  are  in  wavelengths. 

Newton's  method  is  used  to  solve  for  the  intersection  point 
to  an  accuracy  specified  by  the  variable  TOL. 
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SUBROUTINE  TANGIVE (XO , YO , ZO , XA, YA, ZA, BB, H, XI , YI , Z I , XSN, YSN , ZSN) 
PARAMETER ( TOL<* 1 . OE-5 ,MXITER-50) 

No  more  than  MXITER  iterations  will  be  performed 
to  attempt  to  find  t  such  that  g(t)<=TOL 
R»(H**2+BB**2)/2./BB 
RSQ=R**2 
AA-R-BB 

DX-XA-XO 

DY-YA-YO 

DZ-ZA-ZO 

DXSQ=DX**2 

DYSQ=DY**2 

DZSQ»DZ**2 

Calculate  coefficients  of  parametric  equation 
g(t)«  At'>2+Bt+C+  2a*SQRT(Dt*2+Et+F] 

A=DXSQ+DYSQ+DZSQ 
B=2* (XO*DX+YO*DY+ZO*DZ) 

C=XO**2+YO**2+ZO**2-RSQ-AA**2 

D=-DZSQ 

E=-2*Z0*DZ 

F=RSQ-Z0**2 

Solve  for  intersection  using  Newton's  method 


T=.S 

DO  10  1=1, MXITER 
TSQ=T**2 

GC=SQRT ( D*TSQ+E*T+F ) 

GT=A*TSQ+B*T+C+2 . *AA*GG 
GT  is  the  function  g(t) 

GP=2*A*T+B+AA* (2*D*T+E) /GG 

GP  is  its  first  derivative,  g'(t) 

TM=T 

T=TM-GT/GP 

IF(GT.LT.TOL)  GO  TO  20 
CONTINUE 

Calculate  intersection  coordinates 

XI=X0+T* ( XA-XO ) 

YI=YO+T* ( YA-YO) 

ZI=Z0+T* ( ZA-ZO ) 

Some  quantities  become  undefined  if  the  ray  passes 
directly  through  the  tip,  so  to  avoid  that  situation 
a  small  increment  will  be  added  to  XI, YI 
IF(XI.EQ.O. .AND.YI.EQ.O. )  THEN 
XI=XI+.0001 
YI-YI+.0001 
END  IF 

Find  surface  normal  vector  and  normalize  it;  note  that 
it  is  the  INWARD  norward,  hence  the  sign  changes 
XSN=-2*XI 
YSN=-2*YI 

ZSN=-2*ZI * ( 1 . -2  *AA/ SQRT ( RSQ-ZI *  *2 ) ) 
SNM=SQRT(XSN**2+YSN**2+ZSN**2) 

XSN=XSN/SNM 

YSN^YSN/SNM 

ZSN=ZSN/SNM 

RETURN 

END 
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C  SUBROUTINE  RCCONB 

C  This  subroutine  finds  the  point  of  intersection  of  a  ray 
C  with  a  right-circular  cone  radome,  and  the  surface  normal  vector 
C  at  the  point  of  intersection. 

C  The  inputs  are  XO,YO,ZO,  the  coordinates  of  the  ray’s  initial 
C  point,  XA,YA,ZA,  the  coordinates  of  the  ray’s  destination  point 
C  on  the  effective  aperture  outside  the  radome,  and  the  radome  base 
C  (BB)  and  height  (H). 

C  The  point  where  the  ray  intersects  the  radome  is  XI,YI,ZI. 

C  The  components  of  the  surface  normal  vector  are  XSN,YSN,ZSN. 
c  Ail  dimensions  are  in  wavelengths. 

0sCB81|E8asiBEX»GII«»IIISIIIIIEiaaaBHaitll8liailllBIMMSC>MBMBa»l 

SUBROUTINE  RCCONE(XO, YO, ZO, XA, YA, ZA,BB,H,XI, YI,  ZI, XSN, YSN, ZSN) 

BSQ-BB**2 

HSQ=H**2 

DX=XA-XO 

DY-YA-YO 

DZ<=ZA-ZO 

DXSQ=DX*  *2 

DYSQ=DY**2 

DZSQ=DZ**2 

C  Calculate  quadratic  coefficients 

A=DXSQ+DYSQ-DZSQ*BSQ/HSQ 
B=2* (DX*XO+DY*YO+DZ* (H-ZO) *BSQ/HSQ) 

C=XO**2+YO**2- (H-ZO) **2*BSQ/HSQ 

Solve  quadratic  for  intersection  point  Tj  the  solution 
chosen  is  the  one  inside  ZA  >  2  >  ZO 
S=SQRT(B**2-4. *A*C) 

Tl-(“B+S)/2/A 
T2=(-B-S)/2/A 
21=Z0+T1 *DZ 

1F(T1.GE.0. .AND.T1.LE. 1. )  THEN 
T-ti 
ELSE 
T*=T2 
END  IF 

Calculate  intersection  coordinates 
XI=XO+T*DX 
YI=Y0+T*DY 
ZI=ZO+T*DZ 

Some  quantities  become  undefined  if  the  ray  passes 
directly  through  the  tip,  sc  to  avoid  that  situation 
a  small  increment  will  be  added  to  XI, YI 
IF ( XI . EQ . 0 . . AND . YI . EQ . 0 . )  THEN 
XI*XI+.0001 
YI=YI+ . 0001 
END  IF 

Find  surface  normal  vector  and  normalize  it;  the  sign 
is  changed  because  it  is  to  be  the  INWARD  normal 
XSN=-2*XI 
YSN=-2* YI 

ZSN=-2*BSQ*(1.-ZI/H)/H 
C  normalization 

SNM=SQRT ( XSN*  *2  +  YSN**2  +  ZSN*  *2 ) 

XSN=XSN/SNM 

YSN=YSN/SNM 

ZSN=ZSN/SNM 

RETURN 

END 
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ROUTINE  PATTERN  :  DAN  MCGRATH  t  SEP  87  t  MS  FORTRAN 

This  subroutine  calculates  the  radiation  pattern  of  an  array 
antenna  whose  elements  have  arbitrary  locations,  amplitudes  and 
phases. 

Inputs:  X,Y,2  -  Arrays  containing  the  coordinates  of  the  antenna 

elements  (in  wavelenghts) 

N  »  Number  of  antenna  elements 

A  *  Array  containing  element  amplitudes  (voltage) 

P  -  Array  containing  element  phases  (deg) 

PHI  ■  The  spherical  coordinate  phi  angle  for  the 

pattern  cut 

DT  -  The  angle  increment  in  degrees 

TO  -  Initial  theta  angle  in  degrees 

M  -  Number  of  pattern  points  to  be  computed 

EPE  ■  Element  pattern  factor:  antenna  elements  are 

assumed  to  have  a  coeAEPB  factor  (0  for  isotropic) 
NORM  »  Normalization  option:  If  NORM=0  pattern  will  be 
normalized  so  that  its  peak  value  is  0  dB.  If 
NORM=l  pattern  is  normalized  to  N  times  a 
unit  amplitude  (dBi  pattern). 


SUBROUTINE  PATTERN ( X, Y, 2 , A, P, N , TO, DT, M, PHI , FFP, EPE , NORM) 
DIMENSION  X (N) , V (N) , Z (N) , A ( N ) , P (N ) , FFP (M) 

PARAMETER  (PI=3. 141592  359 , DTR' PI/180. , FLOOR=-100) 

FLOOR  is  the  minimum  dB  level  -  calculated  values 
that  are  less  than  FLOOR  will  be  set  to  FLOOR 
SP=SIN ( PHI  * DTR) 

CP=COS (PH I* DTR) 

FKX=0. 

AMX=0. 

WRITE (*,  *) 

Find  Maximum  Element  Amplitude 
DO  10  J»1,N 

IF(ABS(A( J) ) .GT.AMX)  AMX=ABS( A( J) ) 

Calculate  the  pattern  at  M  angles 
DO  100  1=1 ,M 
TH= (TO+DT* ( 1-1 ) ) *DTR 
ST=SIN(TH) 

CT=C0S(TH) 

C  evaluate  element  pattern  factor 

IF (EPE.EQ. 0. )  THEN 
EP=1 . 0 
ELSE 

IF(CT.EQ.0. )  GO  TO  55 
IF(CT.NE.O-)  EP«=ABS(CT)  **EPE 
END  IF 

EXl=ST*CP*2 . *PI 
RX2=ST*SP*2.*PI 
EX3=CT*2 . *PI 
SR=0 .0 
SI=0.0 

C  . for  each  angle,  sum  over  the  array . 

DO  50  J=1,N 

IF(A(J) .EQ.O)  GO  TO  50 
AMP=A( J) *EP 

ARG  =  P(  J) *DTR-X( J)*EX1-Y( J) *EX2-Z{ J) *EX3 
SR=SP.+AMP*COS  ( ARG) 

SI=SI+AMP*SIN(ARG) 

50  CONTINUE 

55  FFP(I)=SR**2+SI**2 

NAL=M-I 

WRITE (*, ' (1H+,5X,I4,A) ' )  NAL, '  POINTS  LEFT* 

100  IF(FFP(I) .GT.FMX)  FMX=FFP(I) 
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C  Convert  to  dB  and  normalize 

IF(NORM. NE.O. AND. NORM. NB. 1)  FMX-1.0 
IF(NORM.EQ.l)  FMX*PLOAT (N) * *2 
IF(FMX.RQ.O. )  FMX-1.0E-10 
DO  150  I-1,M 

C  set  minimum  power  to  prevent  overflow 

IF(FFP(I) .LT.1.0E-10)  FFP ( I ) -1 . 0E-10 
FFP(I)«10.*ALOG10(FFP(I)/FMX) 

IF(FFP(I) .LT. FLOOR)  FFP(I)-FLOOR 
150  CONTINUE 

RETURN 
END 


Function  PHASE  finds  the  angle  argument  of  a 
defined  as  follows:  Real, Imaginary 


complex  number, 
Definition 


0 

<0 

-pi/2 

0 

>0 

pi/2 

<0 

0 

-pi 

>=0 

0 

0 

>0 

1=0 

arctan( im/re) 

<0 

1=0  _ 

arctan(im/re)+pi 

FUNCTION  PHASE (CC) 

COMPLEX  CC 
PI=3. 141592  359 

CR=REAL (CC ) 

CI-AIMAG(CC) 

IF ( CR. EQ. 0 . )  THEN 
IF ( Cl . EQ. 0 . )  THEN 
PHASE=0. 

ELSE 

IF  (Cl .  GT .  0 .  )  PHASE-PI/2. 
IF(CI.LT.O.)  PHASE*-PI/2. 
END  IF 
ELSE 

ARG=CI/CR 
PHASE=ATAN ( ARG ) 

IF(CR.LT.O.)  THEN 

I F ( PHASE . LE . 0 . )  THEN 
PHASE=PHASE+PI 
ELSE 

PHASE=PHASE-PI 
END  IF 
END  IF 
END  IF 
RETURN 
END 


C=»  SUBROUTINE  CTAYLOR  May  1984  FORTRAN 

C==  by  Dan  McGrath,  Rome  Air  Development  Center  (RADC/EEA)  == 

C==  This  subroutine  calculates  the  coefficients  for  a  Taylor  «=* 

C==  amplitude  distribution  for  a  circular  aperture  array  antenna.  == 
C=-  The  Taylor  distribution  holds  a  number  of  sidelobes  nearest  == 
C==-  the  main  beam  to  a  specified  level,  and  the  remainder  roll  off  =  = 
C==  gradually,  NBAR  is  the  number  of  fixed-level  sidelobes,  and  == 
C==  their  amplitude  relative  to  the  main  beam  peak  is  specified  by  *■  = 
C=*  the  variable  SLLDB,  the  sidelobe  level  in  dB,  e.g.  SLLDB=40  == 
C==  for  first  NBAR  sidelobes  40  dB  lower  in  relative  power  than  == 
C==  the  main  beam.  SLLDB  must  be  a  positive  number!  NBAR  must  *=  = 
C==  be  less  than  10.  == 

C==  The  subroutine  must  be  called  repetedly,  once  for  each  == 

C==  array  element,  R  is  the  element's  radius  normalized  to  the  == 
C==  aperture  radius.  The  calculated  amplitude  weight  is  returned  == 
C==  in  the  variable  GP.  NCALL  is  a  flag  used  to  prevent  repeated  == 
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O*  overhead  calculations!  it  must  be  zero  on  the  first  call,  and  =■ 
C»<>  should  be  nonzero  on  successive  calls.  <== 

C==  This  routine  uses  the  function  subprogram  BESSELJ. 

O-  Reference:  "Microwave  Scanning  Antennas,”  Vol.  I,  ■» 

C=«  by  R.C.  Hansen,  Academic  Press,  1966.  «=« 

QsSSS&BBBE8CteSKSa  =  3BCBaaHBHBKBBCBBaBBaSSCBBBBeSBBBVB3BaBB«X«CaiteeceacB 


SUBROUTINE  CTAYLOR ( SLLDB , R , NBAR, GP , NCALL ) 

DIMENSION  RJ0UM(0s9)  ,UN<0:9)  ,FUM(0i9) 

DOUBLE  PRECISION  BESSELJ 
IF ( NCALL. NB.C)  GO  TO  50 

C  The  following  setup  block  only  needs 

C  to  be  executed  on  the  first  call 

c 

C  UN's  are  the  roots  of  Jl(pi*u) 

DATA ( UN (L ) , L=0, 9 )/0. 0,1. 2196699, 2. 2331306, 3. 2383154, 

+  4.2410628,5.2427643,6.2439216,7.2447598,8.2453948,9.2458927/ 
PI=3 . 141592  359 

ETA=10.0** (SLLDB/20) 

DO  10  1=0,9 
X=PI*UN ( I ) 

10  RJ0UM ( I ) =BESSELJ ( 0 , X ) 

FUM(0)=1.0 

A=ATANH(SQRT ( 1 . -1 . /ETA*  *2 ) )/PI 
SIGMSQ=UN ( NBAR} **2/ (A* *2+ (NBAR- . 5 ) * "2 ) 

DO  30  M=1 , 9 
PRN'JM*1 . 0 
PRDEN=1 . 0 
DO  20  N*1,NBAR-1 

C  The  denominator  product  excludes  the  N=M  term 

IF(N.NE.M)  PRDEN=PRDEN*(1.-UN(M)**2/UN(N) **2) 

PRNUM=PRNUM* (l.-UN(M) * *2/ ( SIGMSQ" ( A* *2+ ( N- . 5)**2)  )  ) 

20  CONTINUE 

FUM (M) — RJ0UM (M ) "PRNUM/PRDEN 
30  CONTINUE 

C  After  the  first  call,  calculations  may  start  here 

50  SUM=0 . 0 

P=R*PI 

DO  60  M=0 , NBAR-1 
SDEN=RJ0UM(M)**2 
X«P*UN(M) 

SNUM=FUM(M) "BESSELJ (0,X) 

SUM=SUM+SNUM/SDEN 
60  CONTINUE 

GP=2 . 0"SUM/PI*  *2 

RETURN 

END 

C==  FUNCTION  BESSELJ  Jun  83-Aug  84  FORTRAN4/5  ■== 

C==  by  Boris  Tomasic  and  Dan  McGrath  == 

C=-  Rome  Air  Development  Center,  RADC/EEA  == 

C==  This  function  subprogram  finds  the  value  of  the  == 

C==  Beeael  function  of  the  first  kind,  JN(X)  of  order  N  == 

C==  and  argument  X.  =  = 

C==  0  <=  N  <=  50  ;  0.0  <=  X  <=  100.0 


DOUBLE  PRECISION  FUNCTION  BESSELJ  (N,X) 

DOUBLE  PRECISION  FACT( 300 ) , TRM1 , TRM2 , S , SUM, CC 

IF  (X.GT.10.0)  GO  TO  100 
IF(X.EQ.O.O.AND.N.EQ.O)  SUM=1.0 
IF(X.EQ.0.0. AND.N.EQ.0)  GO  TO  60 
C  For  small  arguments,  ascending  series  formula  is  usea 


59 


nonnonn  a  nan 


FACT( 1 ) “1 . 0 
DO  10  I-2,N+41 

1C  rACT{  I )  “FACT  ( I  —  1 )  *  FLOAT  ( 1-1 ) 

SUM* (X/2. ) **N/FACT(N+1) 

DO  50  K-1,40 

TRM1  =  { X/2  .  >**(N+K)/FACT(K-H) 

TRM2  = ( X/2 . ) *  *K/FACT ( N+K+l ) 

S=TRM1*TRM2* (-1) **K 
SUM*SUM+S 

50  IF(DABS{S) . LT. 1.0E-10.AND. I . GT. 15)  GO  TO  60 
60  BESSELJ-SUM 
RETURN 

For  arguments  larger  than  10.  a  recursive 
difference  equation  solution  is  used.  Ref. 
Abramowitz  it  Stegun  pg.  385 
100  CONTINUE 

IF(X.LT. 40)  L=INT(16. 1* ( 1 . +X* *2*4 . 6875E-5 ) *SQRT ( X ) +30) 
IF(X.GB.40)  L*INT ( 1 . 5*X+83 . ) 

L1=L+1 
FACT(L1 )=0. 

In  this  block  FACT  does  not  represent  factorials 
FACT(L)=1.0D-60 
DO  130  1=2, L 
II=L-I+1 

130  FACT (II) =2 . * FLOAT (II) /X*FACT { II+l ) -FACT { I 1+2 ) 

SUM=0 . 

DO  140  11=3, L, 2 
140  SUM=SUM+FACT(II) 

CC=1 . / ( FACT ( 1 ) +2 . "SUM) 

BESSEL J=CC*  FACT ( N  + 1 ) 

RETURN 

END 


=  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  -  =  =  c  =  =  :  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

==  FUNCTION  ATANH 

==  This  is  the  hyperbolic  arctangent  function.  It  is  =* 

==  calculated  using  the  identity  == 

==  ATANH ( x )  =  .5  *  ln( ( 1+x ) / ( 1-x ) ) 

==  The  argument  is  reoticted  to  |x|  <  1  == 

FUNCTION  ATANH (X) 

R=(l+X)/(1-X) 

ATANH=.5*ALOG(R) 

RETURN 
END 
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